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Abstract. We present a definition of generating functions of canonical relations, which are 
real functions on symmetric symplectic spaces, discussing some conditions for the presence of 
caustics. We show how the actions compose by a neat geometrical formula and are connected to 
I the hamiltonians via a geometrically simple variational principle which determines the classical 

trajectories, discussing the temporal evolution of such "extended hamiltonians" in terms of 
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I Hamilton-Jacobi-type equations. Simplest spaces are treated explicitly. 

Keywords: Symplectic Geometry, Hamiltonian Systems. 



> 
X 



1. Introduction 

The concepts of momentum and position set together on an equal footing define what is known 
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' today as phase space ||19[, geometrically a difFerentiable manifold M endowed with a symplectic form 

CN ! w (nondegeneracy could be relaxed p^|^5|, but here we keep to the symplectic case). This generalizes 

' Hamilton's formalism, as presented in standard textbooks p| [|lT| [^6| . Usually, however, the phase 

• space is first treated as an euclidean 2n-space, or even a cotangent bundle, with canonical coordinates 
{(Pi^Qi)}- While Hamilton's equation can be generically written down for each real differentiable 

• function on any symplectic manifold, other important and familiar features of the canonical formalism 
. do not always extend in an obvious fashion. Furthermore, although Darboux's theorem allows for a 

local canonical coordinatization of general symplectic manifolds, this is not always the most convenient 
one to use. It is therefore desirable to have a formalism which takes into account the specific geometry 
of the phase space, explicitly, as a further step of the "geometrical revenge" in dynamics p^ . 
X ' Actually, the importance of considering nontrivial geometries can be seen from many different per- 

^ . spectives. Although the law of inertia sets euclidean geometry in a priviledged status, once a system 

is subject to nontrivial constraints its proper phase space geometry ceases to be trivial, generically 
[^|l^](2^]. Similarly, for a system which is invariant under a group of symmetries, it is often possible 
to eliminate redundant degrees of freedom and the reduced phase space is also generally nontrivial 
[j|]|2^]|2^]. Alternatively, not to mention relativistic contexts, one may be interested in the classical 
limit of quantum systems, which are non-euclidean, as spin systems [28|[^|. Equivalently, it is conve- 



nient to have formalisms of classical dynamics intimately connected to some important representations 
used in quantum mechanics. For instance, coherent states |l^||2^ or, more pertinent to this work at 



hand, the "center representation" of operators (cartesian Weyl symbols) |22]|27|. Again, it is desirable 
to have these "semiclassical" formulations in nontrivial geometries, but this requires previously well 
defined classical formalisms. 

Particularly, we need well defined generating functions for finite canonical transformations. When 
M is the cotangent bundle over a configuration space Q, these functions naturally take their values 
on Q X Q. When M is a Kahler manifold ||8|[p!7|], whose metric and symplectic form are related via 
its complex structure, one can define generating functions which are bi-holomorphic on Mc x Mc, 
in which case the trajectories are complex. Despite the utility of the complex and the configuration 
space actions in treating a variety of problems, we can greatly benefit from a formalism leading to 
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actions which are real differentiable functions on nontrivial symplectic manifolds, in relationship to 
the real hamiltonian flows obtained from a real variational principle. 

Here we present such a formalism. This work is a generalization on concepts developed for euclidean 



spaces 1 25] 1 27 1, or for linearized neighborhoods of general symplectic manifolds [41|, dating back to 



the original work of Poincare |3C]. The generating functions are real functions on M itself and not 
on a pair of local lagrangian coordinates. An argument of one such function, m £ M, can be viewed 
as the "center" of the canonical transformation m_ i-^ m+ which the function generates implicitly. 
The corresponding variational principle states that, for an appropriate family of paths : [0,t] ^ M 
whose endpoints are "centered" on m, the action 

h{m {t'),t')dt' + i^^iu > (m, t) 



is stationary for a classical trajectory. Here, {^uj}{'m,t) is the symplectic area between the curve i' 
and the geodesic arc, centered on m, returning from to i^(0). 



The restriction on the full foregoing theory is that M be a symmetric symplectic space [13|[17|[20|. 
This means that M admits of a complete affine connection such that every point m G M is the 
isolated fixed point of an involutive symplectomorphism which coincides with geodesic inversion at 
m. Accordingly, the above formulation of the variational principle is invariant at least with respect to 
general transformations on M preserving the affine connection and the symplectic form. 

The generalization from flat to nonflat symmetric symplectic spaces starts with the notion of double 
phase space, DM p^] g. While in the flat case we define the reflection-translation group, from which 



the flat theory develops (§ 2), the equivalent construction for nonflat spaces is not enough and we must 



use the fact that DM is a symplectic groupoid ||3[|[46|. Then, via a "symmetric exponential map", 
we view all structures of DM in (a subset of) the tangent bundle over M, TM. There we define the 
notion of central groupoid (§ 3) . The properties of the pullback symplectic form allows us to see 
the graphs of canonical transformations as lagrangian submanifolds A of (subsets of) TM and define, 
in § 4, (local) generating functions for finite canonical transformations (more generally, relations), 
which are real functions on M. Such functions generate well defined canonical relations only when A 
is a graph over the zero section of TM. For canonical transformations, there is a further graphical 
condition. These considerations are discussed in § 5. 

The composition of two canonical transformations aj : M — > M , as generated by such "central 
actions" , is treated in § 6. Here, again we find that despite their abstract nature, the rule for composing 
them is very simple. If is the central action for Oj, the composed central action for 02(01) is 

/ai A/„2(m) = Stat 

(mi,m2) 

{faiimi) + /a2("T'2) + A{m , rui , 1712)} 

where A(m, mi,m2) is the symplectic area of the geodesic triangle with given midpoints. On general 
symmetric symplectic spaces, the importance of this function was first realized in the context of star 
product quantization , whose euclidean version has long been well established using the idea of 
centers and chords [|4| 1 27 1. On the other hand, the above rule naturally generalizes the result previously 
obtained on iR^" [^][25|. By iterating these compositions, in § 7, we arrive at the central variational 
principle by taking the limit of an infinite number of infinitesimal canonical transformations (§ 8). 
Finally we discuss the temporal evolution of such "finite-time extensions" of hamiltonians describing, 
in § 9, a central version of the Hamilton-Jacobi equation and, more generally, their time derivative 
with respect to any hamiltonian flow, mixing Hamilton-Jacobi with Poisson brackets. 

While presenting these concepts we often make use of particular spaces to illustrate the text. We 
have focused on the simplest two dimensional cases: the euclidean plane M^, the torus T^, the sphere 
S'^ and the noncompact hyperbolic plane H^. Of course, extending these examples to their cartesian 
products M = M^"', S'^ x • • • x S'^, x T^, etc... is reasonably straightforward and has not been 
considered here. We emphasize, however, that the theory in principle applies to generic symmetric 
symplectic spaces of arbitrary (even) dimensions. 
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2. Central Coordinates on Flat Spaces 

Consider the euclidean plane M^, representing the very simplest phase space of a single degree 
of freedom. is a group, under vector addition, and we can also identify the group product as a 
free transitive action of on itself, the translations: T^x) = x + ^. Further, admits a natural 
involution TZo{x) = —x, which represents a reflection, or rotation by vr, through the origin. Together 
with the identity transformation this gives an action of Z2 on M^. We can form the semidirect product: 
RT := Z2 X M^, which is called the reflection-translation group and can also be seen as a normal 
subgroup of the inhomogeneous symplectic group. Then, TZg = {T^TZqT^^) = = TZqT_2x is the 

element of RT which acts as reflection through the point x: lZ^{x') = 2x — x' . {T^ , TZg} satisfy: 

^|*^|/ — ^(l*-!-!/) ) T^R,^ = ^(^_|.|y2) ' '^sT^ — ^(x-g'/2) ' '^s'^x' — 2^2{x-x') (2-1) 

Now, let (x_,x+) e M"^ X be such that x+ = Tlx{x-) <^ x_ = TZg{x+). Then, x G iR^ is 
called the center of the pair (x_,x+). Actually, x = | (x_ + x+), coinciding, for the euclidean metric, 

with the midpoint of the geodesic arc joining x_ to x+. Identifying ^ = x+ — x_, the transformation 
(x_,x+) I— > (x, ^) is a bijection. 

Finally, we see how the group of translations on can be interpreted as a composition of pairs: 
Let x+ = T^,{xx), X- = r_^„(xA)- Then x+ = T^,T^„{x^) = T^^,_^^„^{x^). But noticing that the 

action x+ = T|*(x_) = x_ can be identified by any of the sets {x_,^}, {x+,^}, {x-_,x+} uniquely, 
we can rewrite the composition of translations as: 

(x_ , xa) (xA, x+) = (x_ , x+) (2.2) 

Although (2.2) was derived from the group product on it can be generalized for spaces which are 
not groups. In other words, we can turn the argument around and identify (2.2) as the fundamental 
algebraic structure on M x M. This introduces the concept of groupoid. 



3. Central Groupoids 

Definition 3.1 Let T,AI be spaces. T is called a groupoid over M, denoted T ^ M, if: 

Gd.O) 3 two maps P- , P+ : F M , called the source and target maps, respectively. 

Gd.l) Let (F X F) D F2 := {{i ,i')\P+{i) = P-{i')]. F2 is called the set of composable elements. 

Then, 3 an associative map "0"; F2 F, (7', 7") — > 7' 07", called groupoid composition, or product, 

satisfying P_{i 7") = P^{i) , P+{i 7") = ^+(7") 

Gd.2) 3 an involution i : F ^ F, called inversion, s.t. V7 G F, its unique inverse 7 = 2(7) satisfies: 

7 (7 7') = i if (7, 7') G r2, (7' 7) 7 = 7' if (Y, 7) e r2 

M can be identified with the set of identities inV, Fe := {7 7, or 7 7 | 7 G F}. Every S Fg 
satisfies P_(7e) = P+(7e), but the converse is not necessarily true. 

Examples: The simplest examples of groupoids are: 
i) F is a group G,Te = {e} = M , T2 = G x G 

a) The pair groupoid F = Af x M, P_ and P+ being the first and second projections. In this case, 
Fg — M is the diagonal, inversion is permutation and the compositon is 

(m_, ttia) (mA, m+) = (m_, m+) (3.1) 

which generalizes (2.2). These two examples are complementary in the sense that they have, respec- 
tively, the minimal and the maximal identity spaces possible. 

We now focus on those groupoids F which are also symplectic manifolds and for which groupoid 
composition respects its symplectic structure ||3|[46|. 

Definition 3.2 A groupoid T ^ M is symplectic if {T,ujy) is a symplectic manifold and : 

LGd) M ~ Fg is a submanifold of T, P± are submersions and i:F^F, 0:F2— >F are smooth. In 



4 



RIOS & OZORIO 



this case T is called a Lie Groupoid . 

SGd) The graph of "0" is a lagrangian submanifold ofTxTxT, where T ^ T is an antisymplectic 
isomorphism: (f,u;p) = (F,— wr)- In shorthand notation, we write this as an "additive" property : 

a;r(7i0 72) ~ wr(7i) + '^r(72) • 

As consequences, M ~ Fg is a lagrangian submanifold of T and inversion is antisymplectic, i.e. 
the graph of "i" is a lagrangian submanifold o/F x F , wr(7) ^ — i^r(7) • 

Example: let M be a sympletic manifold, uj its symplectic structure. Then, the pair groupoid M x M, 
with sympletic structure 1^x^ = 810, where 

5uj:=PI{uj)-P1{uj)=uj+Quj^ , (3.2) 

is a symplectic groupoid, called the product or double phase space and denoted DM. That is DM = 
(M X M) ^ M. Notice that while DM satisfies SGd, the cartesian product M x M, with the usual 
symplectic form w © u, does not. 

In order to define the notion of "center" precisely, we need to consider phase spaces which are 



symplectic, but also have the property of being symmetric |13| |17| |20|: 



Definition 3.3 A differentiable manifold M endowed with a symplectic form uj shall be called a 
symmetric symplectic space if (M, lo) admits of a complete affine connection V which is symplectic, 
ie. Vw = 0, such that Vm G M, m is the isolated fixed point of an involutive symplectomorphism 
TZm : M ^ M which coincides with the inversion, at m, of the geodesic flow, for all geodesies through 
m. Accordingly, TZm satisfies TZmT^m'T^m = ^7^™(m,') ; {T^mT^m'} is the group of displacements on M 
and V is the corresponding canonical connection with null torsion and parallel curvature. 

This generalizes the refiection-translation structure on euclidean space. An important particular 
case is when V is the Levi-Civitta connection associated to a complete riemannian metric rj. In 



this case, TZ^{r],uj) = {r]^uj) and one can show |15|[34| that {M,uj,ri) is a hermitian symmetric 



space, that is, M is a Kahler manifold Qll?] whose complex structure J satisfies VJ = and 
r]{JX,Y) = uj{X,Y), yX,Y £ X{M) , and for which the involution TZm is holomorphic. More 
generally, part of the following theory, namely the local definition of central generating functions (§ 



4), can be developed on general symplectic manifolds |3C]|4C]|41]|42]. However, the geometrical law 
of composition (§ 6) and its consequences need the full setting of a symmetric symplectic space. As 
simplest examples of these, we shall consider the flat euclidean plane and torus T^, both being 
also groups, and as simplest nonflat examples, the sphere 5^ and the noncompact hyperbolic plane 
H^, both being hermitian symmetric spaces. 

When M is a symmetric symplectic space, let us denote the exponential map by Expm : T^M M , 
V ^ /Or(l) 1 where Prit) is the geodesic in M defined by Pr(0) = m , PriO) = v . 
We can define a (restricted) groupoid structure directly on (a subset of) its tangent bundle TM, 
associated to DM via the symmetric exponential map : 

^ : TM^DM , r = {m,v) ^ 7 = {Exp^{-v) , Exp„(t;)) . 

Here, <I> is "symmetric" in the sense that (7^mOcI))(rn,, v) = ^(m, —v) and thus 7^7 = ^(7). Denoting 
Q = ^*(6uj), 6(jj defined in (3.2), we have that is a closed 2-form on TM since 5u) is a symplectic form 
on DM. If M has no closed geodesies, as for a hermitian symmetric space of noncompact type (or M^"'), 
then $ is a bijection and (TM, Q) is a sympletic manifold. But otherwise, $ is not globally invertible. 
Consider then the set (TM)o C TM, {TM)q being the maximal subset, connected to the zero section 
r°M = {(m,0)} C TM, on which $ is invertible. Denote ^>o := ^|(TAf)o' ^0 •= ^l(TAf)o = ^o('^'^)- 
Then ((TM)o, Qq) is a symplectic manifold. Now, it is known that for complete riemannian manifolds, 
only the stratum of Si C TM has maximal dimension, where ai{m) = n T^M is the subset of 
TmM defined by I{v) = i , i = 1,2, ... , with I(v) = number of geodesies of length ||f || which connect 
m and Expm{v) , see Thus, if M is a hermitian symmetric space of compact type (or the torus), 
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although $0 : iTM)o — > DM is not onto, it only leaves out a set of measure zero with respect to the 
Liouville form on DM obtained from Su. We conjecture that this may be true for other symmetric 
symplcctic spaces with closed geodesies, as well. These facts motivate a slight modification in the 
concept of symplectic groupoids, suitable to our purposes: 

Definition 3.4 Let M he a symmetric symplectic space. {(TM)q,Qq) shall be called the standard 
central groupoid over M, denoted {TM)q ^ M , satisfying : 

CG.O) 3 three maps Po,P-,P+ : {TM)q — M, called the central, source and target maps, respec- 
tively, where Po{m,v) = m is the natural projection and P±{m,v) = Exp^{±v) . 
CG.l) On the restricted set of composable elements ((TM)o x (TM)o) D (TM)^ := {(r',r") | 
P+(r') = P-(r") and (P-(r'), P+ir")) £ Image{^Q) C DM}, the composition "0"; {TM)l 
{TM)o, satisfying (Gd.l), is defined by r' r" = (P_(r'), P+(r")) • 

CG.2) 3 an involution i : {TM)q {TM)q, called inversion, satisfying all properties in (Gd.2) plus 
Pq {i{T)) = Po{t), Vr G {TM)q. Specifically, if t = {m,v) then f = i(r) = (m, —H) . 
CG.3) The identiy space is the zero section and inversion is smooth. However, P± are only locally 
submersions and is locally smooth. {TM)o is a local Lie groupoid . 

CG.4) The graph of "0" is the union of lagrangian suhmanifolds of {TM)q x (TA4)q x {TM)q. The 
graph of "i" is a lagrangian submanifold of {TM)q x {TM)o, the zero section being lagrangian in 
{TM)o. We may call ((rM)o,rio) o, local symplectic groupoid . 

When M has no closed geodesies, {TM)o = TM, $o = TM is diffeomorphic to DM. In this 
case, the words "restricted" and "local" do not apply and the central groupoid is a bonafide symplectic 
groupoid with the extra structure of a central map. 

Examples: {TlR^)o = TM^, (TH\ = TH^, the trivial cases since both and H"^ have no closed 
geodesies. For we have the following: Let AS^ denote the set of antipodals in x S^, i.e. 
{((/?, ± tt)}, under the usual representation 3 e*^, if G [0, 27r]. The subset TT^ D (T7"^)o := 
{r= {p,q ; Vp,Vq)\ \vp\,\vq\ < 7r/2} is s.t. $o : {TT'^)o DT\{AS'^ x AS'^) is bijective. Clearly 
AS^ X AS^ is a set of measure zero in T'^ xT'^. For 5^ the situation is very similar: Let AS"^ denote 
the set of antipodals in DS^. Taking the subset TS"^ D (TS'^)o := {r = {m,v)\ \v\ < 7r/2} , then 
$0 • iTS'^)o DS'^\AS'^ is bijective. Again, AS'^ has measure zero in DS^. 

However, when M has closed geodesies, the standard central groupoid is not the only central 
groupoid possible. To see that, consider another subset {TM)i C TM which satisfies: {TM)i n 
(TM)o = and $((rM)i) = $((TM)o) = iDM)o C DM . Furthermore, denoting the map 
restriction by $i = we limit attention to those subsets on which $i is almost everywhere 

injective. In other words, the set A^C {DM)q on which is multiple valued has measure zero with 
respect to the Liouville form on DM obtained from 6ui. Then, defining an equivalence relation among 
different pre-images in (TM)i of the same point in {DM)o, we have another bijection $i : (TM)i — 
(DM)o, where {TM)i = {TM)i/ ~i with the equivalence relation r ~i r' iff $i(r) = $i(t'). 

Accordingly, we denote such points in (TM)i by f = [r]i, where [r]i = {r' G {TM)\ \ r' ~i r}. If 
M C {DM)o is the diagonal, its pre-image $^^(M) =: T^M is a global cross section of (TM)i setting 

a fibration Pi : {TM)i THL s.t., if f G (fM)i and G THd then Pi(f) = iff Pq{t) = Po(r^). 
In other words, the central map Pi is essentialy Pq : TM — > M itself. We denote Oi = ^\{5uj). 

Definition 3.5 We call {{TM)i,^i) {THL,uj) a nonstandard central groupoid. On {(fM)iX 
{fM)i) D {fM)l := {(f,f')| P+{f) = P-{r') and^{P^{f),J>+{T')) G {DM)o}, where P±(f) = 
P±iT)} G ^; th,e restricted groupoid product : {TM)1 (TM)i is given by fQf' = ^i^{^i{f) 
4i(f')), and inversion by i{f) = 4j^^(i(^>i(f))). // {TM)2 T^M is another central groupoid, it 
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is strongly equivalent to (TM)i if, for every f e {TM)i and f G {TM)2 s.t. I>i(f) = ^2{f'), it 
follows Pi(f) = P2{f'). //$i(f) = 1*2(7") but Pi{f) P2{t'), f' is weakly equivalent to f. 

Examples: When M = T^, there are infinite central groupoids, but only 4 strongly inequivalent ones: 
(rT^)oi = {r = {p,q,Vp,Vq)\ \vp\ < tt/2 , \vq — 7r| < 7r/2}, similarly for (TT^)io, where \vp — 7r| < 
7r/2 , \vq\ < 7t/2, and (TT^)ii, where |fp — vrj, — 7r| < 7r/2 , together with the standard one (TT^)o = 
(T7"^)oo- In this case, all the nonstandard central groupoids are simple subsets of TM, i.e. no 
equivalence relations had to be considered since each = ^\(^rpq-2^.. is already injective. Each point in 
(Z)T^)o has 4 strongly inequivalent pre-images: (j),q,Vq,Vq), {p±7T ,q, Vp^w, Vq) , (p, Q'lbTr, VpjVq^ir) 
and (p lb TT , g ± TT , Vp ± n , Vq ± n). In the case of 5^, there are only 2 strongly inequivalent 

central groupoids. The standard one and a nonstandard central groupoid: (TS'^)i = {TS'^)i/ ~i , 
where {TS'^)i = {{m,v)\\v\ € ('7r/2 ,vr]}. The equivalence relation is nontrivial only when \v\ = vr, 

coinciding with the definition of the nonstandard pre-image of the diagonal, T^S"^ := { [r^]i , where 
= {m,vi) G {TS^)i is s.t \vi\ = it and {m,Vi) ~i {m',v[) iff m = m' }. Thus, each point on 
(-DS'^)o has two strongly inequivalent prc-imagcs. If wc denote ^Q^{m^,m^) = {m,v) = {9,(p;v, (3), 
see below for definition of the fiber coordinates {v,(3), then i^^(m_,m+) = (tt — 0, — 7r;7r — — tt). 

The reason for considering such nonstandard pre-images in TAI is, of course, that when M has 
closed geodesies, the geodesic arc connecting two points in M is not unique. Accordingly, the standard 
central groupoid refers to the shortest of such geodesic arcs and the /^-nonstandard central groupoid 
refers to the fc-fold ones. Here, a fc-fold geodesic from m' to m" is defined as the geodesic (m' — >■ 
m) * Lk{m) * (m rn"), where Lj^{m) is a fc-fold geodesic loop based at m, the midpoint of the short 
geodesic, and * denotes free product. (If k = 0, Lq = id is the trivial loop and {m' m") is the short 
geodesic). For instances of nontrivial geodesic loops, on T^, = Lf^g^ ^ = * L^, where Lp, Lq are 
the single irreducible circuits, while on S'^, Lj^ is a |fe|-repetition of a geodesic meridian and it is easy 
to see that in this case we only need to consider strongly inequivalent geodesies. We also consider : 

Definition 3.6 Let {TM)f- := {(f,f') G {fM)i x {fM)j \ P+(f) = P_(f') and (P_(f), P+(f')) 
G {DM)q}. We define generalized, or mixed compositions (which are still restricted, in the original 

sense) as c^j : {TM)fj {TM)k , by cf^(f,f') = i^^($j(f) ij(f')). If M is not simply connected, 
we may want to consider only mixed compositions whose geodesic triangles of composition are reducible 
circuits, in which case the triple {i,j] k) is called reducible. When M is simply connected, it is enough 
to consider mixed compositions within the set of strongly inequivalent central groupoids. 

Now, the application of central groupoids in classical dynamics rests on the following : 

Definition 3.7 Let N be a symplectic space. A polarization on N is an intergrable foliation of N 

by lagrangian leaves . 

Lemma 3.1 Pq sets a polarization on {(TM)q,Qq) , called the vertical or central polarization . 

Proof: First, every fibration is an integrable foliation. Second, since each fiber is (an open subset 
of) the tangent space at each point in M, the fiber dimension is half the total dimension of (TM)o. 
To prove isotropy with respect to = ^o('^'^)' iiote that inversion is anti-symplectic in (TM)q, 
i.e. i*{^o) = —^0- But, denoting the fiber restriction by (TmM)o, we have that i^ := j^^y^ = 

$Q ^ o TZm ° ^0- Since uj is TZm invariant, i.e. 'R.*^{uj) = lo, it follows that i^ ((^^o)m) = (^o)m) where 

(f^o)m = ^o\(T^M)o- (f^o)m = . □ 

Further : we realize that Lemma 3.1 applies to every Pi : (TM)i — ^ T^M, as well . 

Thus, in what follows we should consider all possible central groupoids similarly. However, we shall 
mostly deal with the standard central groupoid, both to simplify the treatment and because we'll 
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often focus on the limit 1^1 0, which can only take place in {TM)q. Accordingly, we often drop the 
denomination "standard" in following definitions and discussions. 

Finally, we should compare the central groupoids with the local isomorphism DM T*M, in a 
neighborhood of the diagonal in DM, for any symplectic manifold M [^] [^] The differences lie in 
the scope and properties of the symplectic structure. While {{TM)i,Qi) contains explicit information 
on the affine geometry (geodesic structure) of M, the canonical symplectic form on T*M does not. 

Formulae: To base the abstract concepts of this paper, we'll often illustrate them in the simplest 
spaces M = , , . We now provide some useful local formulae for their standard central 

groupoids. Start with the flat spaces, locally identical. Take coordinates on as canonical pairs, 
with usual metric and symplectic form. The tangent bundle is also a linear space : TM^ 3 f = 
(x ; v) = (p, q ; Vp, Vq) , and the symmetric exponential map is : 

<l>(r) = <I>(x ; v) = {x — V ; x + v) = (x_, x+) (3.3) 

with pull-back symplectic form given by 

= 2dTj , S = Vpdq — Vgdp . (3-4) 

For the sphere, take local polar coordinates: S'^ B m = {6,ip) , ^ G [0, tt] , € [0, 2it] . We shall be 
using the following abbreviations: 

Sa = sin{a) , Ca = cos{a) , Ta = tan{a), with S~^{f),C~^{f),T~^{f) 

denoting their respective inverses. With usual metric and symplectic form, the natural coordinates on 
the tangent bundle are TS'^ 3T = {e,ip; e,ip), but it is more convenient to introduce polar coordinates 

on the fibers as well, t = {6,(p ; v, 0) where: v = \v\ = \J 0"^ + Sg^p^ , vCp = 9 , vSjs = Sgip . With 
these local coordinates we write the symmetric exponential map as 

Lp± = ip± T-'^{SySp/g±) + (1 - Sign{g±))TT/2 , g± = SeCy ± CgSyC/s 

and the pull-back symplectic form 0,q = ^q{6uj) is given by 

J7o = 2dS , S = Sy{CpSed^ - SpdO) . (3.6) 

On , we adapt the local spherical formulas by letting 

I — > ip , V I — > ip (3.7) 

where i = \/— T. Thus, i — > iSp, Cq i — > Cp , 5„ i — > iS^, i — > Cp , with 

Sa = sinh{a) , Ca = cosh{a) ,fa = tanh{a), and S^^{f), C^^{f), f~^{f) 

the respective inverses. Then, from (3.5) and (3.6) we get <I> and —Q, for H^. 

We can see explicitly from (3.4), (3.6) and (3.7) that the vertical spaces are isotropic in these 
examples. The same holding for the zero section T^M ~ M. Notice also that the pull-back symplectic 
form is the exact derivative of a symplectic potential without any vertical differential components. 
This fact shall be thoroughly exploited in what follows. 

4. The Central Equation 

We started by emphasizing the algebraic structure on double phase spaces, or on central groupoids. 
Historically, the symplectic structure was predominant, however, for it introduced the very useful 
concept of action, or generating function of a canonical transformation on the original, or simple 
phase space {M,lo). Let a be such a symplectomorphism M M, a*(t(;) = w. Its graph Ca in the 
double phase space, DM D Ca ■= {(m_,m+) | m+ = a(m_)}, is a lagrangian submanifold for the 
symplectic form 5u;, i.e. 5Lo\Ca = 0. Similarly, if Ca C Image(<l>o) = {DM)o C DM, then its pre- 
image := $q ^(£q) is lagrangian in the central groupoid {{TM)q, Qq). Conversely, every lagrangian 
submanifold A^ := ^Q^{Ca) defines a symplectic or canonical relation on M [^|^4|, which is a 
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canonical transformation when is a graph over M_ C DM. Genericahy, we consider those subsets 

C £ai satisfying C Image(#o)> and their corresponding pre-images A^^ in (rM)o. 
Now, the definition of local generating functions for lagrangian submanifolds depends on the choice 
of a local symplectic potential (a local 1-form whose derivative is the symplectic form) suitable to a 
polarization which is, at least locally, a fibration over a referential lagrangian submanifold containing 
the supports of those generating functions [37|[^][^][|^][^]. For the central groupoid over M, we 



take M itself, also seen as the zero section, as the referential lagrangian submanifold corresponding to 
the central polarization. The suitable symplectic potential is given by : 

Definition 4.1 Let {(TM)q,0,q) be the standard central groupoid over (M,uj). A symplectic potential 
Zq, for r^O; shall be called a standard central potential if it satisfies 

X\Zo = P!i{X)\Zo , yX G X{{TM)o) . (4.1) 

Here, Pq is the differential of Pq and x\ a denotes the vector-form contraction. 

Of course, any Z = Zq + dQ, Q G C|j((TM)o), is another potential, but not generally central. 
Condition (4.1) tells us that central potentials have no vertical differential components, so we can 
identify these potentials explicitly in our examples as Zq = 2S, from equations (3.4), (3.6) and (3.7). 

In these particular examples, Zq is a global potential on {TM)q i.e., VIq = dZQ is exact. We argue 
that the general case follow these known examples. First, notice that such a potential always exists 
in a small neighborhood of the zero section T^M. For, take the exact 2-form uj on TM defined 
by d}{m,v) := d{v\uj), ViT G X{M). Then, uj = L^{w), the Lie derivative of w, since uj is closed. 
On the other hand, when |u| = t ^ , QQ{m.,v) ^ 21im ^ /q L^(u;)(it'} ~ 2L^{lo). That is 

Oq — > 2u; = 2d{v\uj). Thus, in a neighborhood of T^M, JIq ~^ where C,{m,v) := v\Li!{m) clearly 
satisfies (4.1). (on the flat examples, Zq = 2C, globally, but generally 2^ only approximates the central 
potential in a small neighborhood of T^M ). Finally, T^M is lagrangian for Q.q and each vertical fiber 
is also lagrangian, as well as contractible. Thus, VLq = dZQ is exact. 

Recalling that ^Zo(m,iT) ^ ({m,v) = v\oj, as l^l 0, we may refer to the 1-form ^ as Hamilton's 
potential, since it is intimately connected to Hamilton's equation. To see this, remember that we can 
write the latter as a map z[h] : M — > TM, V/i G C^{M), by m i-^ Th, where Th = {m,Vh) and Vh 
is given by {dh + Vh\uj){m) = 0, or equivalently, C('^/i) = —dh{m), Vm G M. Or still, by denoting 
TM D Afi := graph of z[h], hamilton's equation becomes Clj^h ~ implicitly defining z[h]. 

Similarly, the existence of a central potential for il.Q allows for a "finite time" extension of Hamilton's 
formalism , now in the context of generating functions. Thus, suppose that a lagrangian submanifold 
Aq, C (TM)o is locally a graph over r°M ~ M, that is: 

Rank(P^ | r^A„) = 2n = dim{M) , Vr G A^, s.t. Pq{t) = m G U C M (4.2) 

then, since dZQ\Aa = 0, from (4.1) we obtain: 

Proposition 4.1 For every lagrangian submanifold Aq, C {TM)q satisfying (4-2) there exists a stan- 
dard central generating function fa G C|j(M), satisfying the standard central equation 

= dfa , on {TU)q C {TM)q . (4.3) 

Ac 

Conversely, for a given fa, the above equation defines Aq, implicitly, i.e. it provides a standard section 
z[fa] = Fa : U — > (TU)q which is well defined when condition (4-2) is satisfied and as long as fa 
satisfies appropriate standard consistency conditions . 

We shall see examples of such consistency conditions shortly. Now, the map Fa can be multiple 
valued if {TmM)Q n A Q is not unique, in which case we should break Aq, into branches Aq'', each one 
uniquely given by a map Fa '^ in a subset of M via Zq\ (r) = dfa \ for each fa \ except for {t^} C Aq 
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where Rank (Pq | Tr^.f^a) < 2n. Then, via the symmetric exponential map, (4-3) locally generates a 
lagrangian submanifold Ca C DM, a canonical relation on M. 

The similarity between Hamilton's equation and the central equation is striking . However, the former 
provides infinitesimal transformations while the latter generates finite relations. For this reason, 
not every canonical relation can be generated by a real function on M via the central equation, 
everywhere. Generically, the presence of central catastrophes, {rk} C Aq s.t. Rank (Pq | Tr^Aa) < 2n, 
is unavoidable. Their projections, {m^} = {PoiTk)} C M, are called central caustics. To circumvent 



this problem, new sets of "complementary" generating functions are needed |25|. Alternatively, we 
can let the functions depend on extra parameters and look for their stationary points 

On the other hand, by introducing a real parameter X = t/2 (a scale) in the map z[h], for instance, 
by multiplying every hamiltonian h by t/2, we can see Hamilton's equation as a map from M into a 
small neighborhood of the zero section in TM, if t is sufficiently small. Since in such a neighborhood 
^Zq and ( are approximately equal, fa = —th is a central generating function for the infinitesimal 
canonical transformation generated by h via Hamilton's equation. Since these are always well defined, 
every infinitesimal canonical transformation can be generated by a central function. This is obvious 
if we notice that such transformations are small deformations of the identity, associated to lagrangian 
submanifolds in (TM)q which are small deformations of the zero section, thus satisfying (4.2). 

Furthermore, the map = —tH : M {TM)q , obtained via the central equation from function 
/q, = —th, where h is the hamiltonian, is given by m i-^ {m, ^Vh{m)), for sufficiently short times t ~ e. 
Since ^Vh — > , as e — > , for very short times t ~ e — > 0, the implicit canonical transformation 
(m_,m+) obtained via the symmetric exponential map coincides with the linearized version (x — 
^Vh,x + ^Vh), for any choice of local linear coordinates on [/ C M, s.t. m ~ x, regardless of the 
specific affine geometry of M. Thus, for very short time motion, the geodesic segment centered on m 
converges onto the hamiltonian orbit that propagates from m forwards and backwards in time. 

Now, the pertinent remark should be made that the concept of generating function is traditionally 
defined on the double phase space DM, but this usually requires us to previously take polarizations on 
M itself. For general cotangent bundles, with vertical polarization, the generating functions naturally 
take their values from pairs of base space points, / = f{q-,q+)- For general Kahler manifolds, with 
complex polarizations, the natural generating functions are bi-holomorphic functions / = f{z-, 2;+), z± 
complex. The present approach allows us to consider real generating functions on a general symmetric 
symplectic space M itself, defined via a real polarization on the central groupoid. These functions 
can sometimes be pictured as 'finite time hamiltonians" . More generally, such central generating 
functions can locally be defined on any symplectic manifold by considering a local polarization of DM 
which is transversal to the diagonal M C DM. If linear coordinates are chosen on a neighborhood of 
a point in M, these functions can be mapped to Poincare's generating functions However, the 

full geometrical properties of the central generating functions need the symmetric symplectic setting. 

Also, we must point out that other versions of the central equation are available : 

Proposition 4.1' When M has closed geodesies, we can choose nonstandard central potentials Zi s.t. 

dZi = ^*{6uj) and X\Zi = Pl{X)\Zi , e X{{TM),) . (4.1') 
Each lagrangian submanifold = ^i^{J~-a) C {TM)i, satisfying 

Rank{Pl\TrAi) = 2n , over U C M TH/I , (4.2') 
is generated by a nonstandard central function f^ via a nonstandard central equation : 

Z^\^^= dfi , on {fu)i C (fM)i , (4.3') 

which provides a nonstandard section : U ^ {TU)i defining A^. When no distinction is made, or 
by means of generalization, we refer to (4-3) or (4-3') simply as the central equation . 
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Notice that we have abbreviated the standard notation: = fa, A.^ = Aa, = Fa, in previous 
definitions. Also, the index "z" in (4.3') is not the same as the index (r) defined earlier, referring 

to difi'erent branches of a single pre-image of the set {c'^a^} C n (L>M)o. Thus, a generic central 
generating function can carry up to two indices {i, (r)}, in order to be fully identified. 

Finally, since central potentials are defined modulo exact differentials on M, it is important to 
emphasize that each nonstandard central equation, with its appropriate nonstandard consistency con- 
ditions, is defined with respect to a choice of the corresponding nonstandard central potential. We may 
fix the standard one by setting ZqIj^om = , but similar choices for the nonstandard ones are not 
necessarily the best. In ^ 6 , Proposition 6.1' , we describe a consistent choice for all Zi's which is 
suitable for the mixed composition of central generating functions . Such compositions, as we shall see, 
exhibit very neatly the full symplectic and affine geometry of M itself, a fact having some interesting 
bearings on the problems of quantization and semiclassical analysis. 

Examples: Standard cases only. We show the explicit map z[fa] = Fa : M ^ {TM)q, given by each 
generating function /„, and the canonical relation (m_,m+) G DM. To get rid of factors of 2, we 
often rescale and identify fa = 2/. On IR^, using (3.4) and (4.3), F is written as 

Vp = df/dq , V, = -df/dp , or ia=2v = -J- [dfa/dx] . (4.4) 

where J is the symplectic matrix on M^. Composing with (3.3), we have the canonical relation as : 

x± = x=fJ ■ [df/dx] = X =F ^ J • [dfa/dx] (4.5) 



See |25||27|. On the torus, (4.4)-(4.5) are valid, but we impose a standard consistency condition which 
ammounts to constraining the map (4.4) onto (TT^)o only : 

\df/dq[ , [df/dp[ < 7r/2 , on [/ C . (4.6) 

On S^, by (3.6) and (4.3), F : S"^ — > (TS'^)o is written in polar coordinates as : 

v = S-Hs'{f)) , p = -T-'S^^^^^+{l-Sign{df/d^))7T/2 (4.7) 

where, using the contravariant metric on S'^, we define the symbols 



s\f) := ^ ^{df/dey + {{i/Se){df/d^)y , cO(/) := Vi - {s'ifW ■ 

Notice that (4.7) has real solutions only if / satisfies the consistency condition 

S\f)^[[df[[<l , onUcS'. (4.8) 
Composing with the symmetric exponential map (3.5), we obtain 



0± = C-'{CeC^{f)Tdf/dip) \ 
^± = VT T'^Se{df/de)/a±} + (1 - Sign{a±))7r/2 , / 



(4.9) 



where a± = SgC^{f) ±Cg{df /dip) , as the local expression for the canonical relation (m-_,m-|_) which 

-•k ( 
•-iRV 



is generated by fa = '^f ^ ^Mti^"^) i satisfying condition (4.8). Similarly on H^, by (3.6), (3.7) and 



the central equation, the map F : H TH is given in local polar coordinates as : 

^ = S-\s\f)) , /3 = -r-^ I ^f^^'g^^ I + (1 - Sign{df/d^))n/2 (4.10) 
which now requires no consistency condition. Once more, we have defined 



s\f) ■■= \\df[[ ^ ^J{df/dpY + {{i/~s,){df/d^)Y , c\f) := 0T(^W 
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Finally, composing with the symmetric exponential map we have the local expression 
P± = C-\Cp&{f)±df/dip) 

V± = ^^ T-i [Sp[df/dp)h±] + (1 - 5i5n(7±))7r/2 , 
where 7± = SpC^{f) ± Cp{df/d(p) , for the corresponding canonical relation on H^. 

5. Central Actions and Relations 

We saw in the last paragraph that a function G C^{M) can locally be taken as central generating 
function of a canonical relation C {TM)i, provided (4.2') holds. This is a "graphical" condition over 

U d M T^M c {TM)i and is therefore written with respect to A^. But generically, it is precisely 
this submanifold that needs to be found given /q. Furthermore, it is important to distinguish which 
of the various lagrangian submanifolds in (TM)i do correspond to canonical transformations on M, 
i.e. which are pre-images, under the symmetric exponential map, of graphs over M_ C DM. Again, 
since we usually start with the central generating function, from a practical point of view we need 
such a distinction, as well as an alternative to condition (4.2'), written directly in terms of We 
now proceed in this direction, in the standard case. First we obtain : 

Lemma 5.1 A function f = ^fa G C^(M), k > 2, can locally he the standard central generating 
function of a canonical relation, via central equation, only if it satisfies all consistency conditions 
required for the definition of the map F : M D U ^ {TU)q, F : rm-^ v E {T^M)q, and 

\det[dFydm^]\ < oo , Vm e C M , (5.1) 

for any choice of local coordinates {m*} on U, {v*} on (r^M)o, with {F*(m) = v'^{t)}. 

Proof: Let Tq G Aq, C {TM)q, Po(7'a) = ^ ^'^d take local coordinates {a;^,-- - ,x^^} on 
a neighborhood X C A^ of the point r^, {m},--- ,m^"'} on a neighborhood U <Z M oi the point 
nia and {v^,--- on a neighborhood V C (T^^M)o of the point = Pv{Ta), P-u '■ N ^ V , 

N = {U xV) C {TM)o . Then, any point t e X C N is locally written as r = {x*(r)} = {x'}, 
T = {m^{Po{T)), t;*(P„(T))}, or r = {m^ {Pq{{x^)) , v^{Py{{x^}))}. Now, the graphical condition 
(4.2) is simply [dm*] = Aq ■ [dx^ , [^q ] = [dm^dx^ , satisfying < \det{Ao)\ < oo , on r, that 
is [dx*] = • [dm^] , < \det{AQ^)\ < oo , and the failure of (4.2) is written as \det{Ac)\ = 
<^=4> \det{AQ^)\ = oo . But since no graphical condition over V is assumed, we have only [dv'^] = 
Ay ■ [dx^] , [Ai!'v] = [dv'^ /dx^] , satisfying |def(A„)| < oo , on r, and the same cannot be said of A~^ 
because |fiet(ylt,,)| = is a real possibility. Thus, \det{A.^ ■ Aq^)\ = \det[dv'^ /dm^]\ < oo is a necessary 
condition for the existence of central generating functions. Now, if such exist, then they satisfy (4.3) 
and we can write the map F generated by / G C^{U) as {v'^{r) = F*(m)}, provided / satisfies any 
required central consistency condition. □ 

Examples: On M^, (5.1) becomes |(iet[5^/]| = |(iet[9^//5x*5x-']| < oo . On we get the condition 
\det[d'^f] + {Ce/Se)idf/dip){d'^f/dedip)\ <oo ,onU CS^ , where det[d'^f] = (5V/56'^)(5V/ V) - 
(d"^ f /dOdif)"^ , provided {6 = 0,7r) ^ U, and / satisfies < \\df\\ < 1. The restriction on 6 is 
easily removable by choosing a new origin for the polar coordinates on S"^, but a new local anal- 
ysis is needed when ||d/||(mo) = 0. Since — implies \v\ — ^ 0, we know these points cor- 
repond to fixed points of any canonical transformation which can be generated by /. Thus, a 
(local) alternative consists of expanding / around ijiq in linear coordinates and applying the pre- 
vious flat equation, on a small neighborhood of mo- Similarly, we can use that ~ v to get 
\det[d'^f] + {Cg/S0){{df/d(p){d'^f/ded^)-{df/de){d'^f/d(p'^)}\ < oo , as the local form of (5.1) 
on a smaller neighborhood U' of a point niQ £ U for which ||(i/||(mo) = . Similarly, on [7 C H^, 
(p = 0) ^ U, we write (5.1) locally as \det[d'^ f] + {Cp/ Sp){df /d^){d'^ f /dpdf)\-{SpS°{f)C°{f))-'^ < oo, 
if > 0, and \det[d^f] + {Cp/S,) {{df/d^){d^f/dpdif) - {df /dp){d^ f /d^^)] I • (Sp)-2 < oo on 
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a smaller neighborhood of the point mo for which ||(i/||(mo) = 0. Or we can expand / around mo in 
linear coordinates and use the fiat equation. 

Lemma 5.1 provides a necessary, but not sufficient condition for the existence of central generating 
functions, generically. In other words, when A C {TM)q is a generic canonical relation, it is possible 
that |(fet(j4o)| = |(iet(^t,)| = 0, on r € X C A. That is, we don't necessarily have a graph over 
either M or Tj^M, on r = (m, v) G A, generically. Even when restricting to canonical transformations 
on M, it is possible that det{Ao) = det{Av) = 0, if dim{M) > 4. To see this, consider the simple 
example: M = IF^ x iR^, u = U{i^ © a;(2), a = "T^q^' ® id^'^\ Furthemore, (5.1) does not distinguish 
between canonical relations and transformations. In this respect, we extend Lemma 5.1 to 

Lemma 5.2 A function f = ^fa G C^{M), k > 2, can locally be the standard central generating 
function of a canonical transformation on M only if it satisfies condition (5.1), besides any central 
consistency condition required for the definition of the map F : M D U ^ {TU)q, m ^ F{m) G 
(T^M)o, via the central equation, and 

< \det[d{Expm{-F{rn)y /dm^W < oo , e U C M , (5.2) 
for any choice of local coordinates on M . 

Proof: Consider Aq, = $^^(£q), where £q, is a lagrangian graph over M_ CDM. Let 7 = $o(r) G y C 
La, Y a neighborhood of 7^ = $o(tq;) and consider local coordinates {y^, ■ ■ ■ y^"'} on Y. Similarly, take 
{mi, • • • m^"} as local coordinates on a neighborhood U- C M_ of the point = P-{Ta) = P_(7q). 
Then jC^ is locally a graph over M_ iff [dml_] = B^[dy^], where [B^l] = [dm^_/dy^ satisfies < 
\det{B^)\ < 00 . Otherwise, if £a is not a graph over M_, then from the definition of S_ , det{B-) = 0. 
But since $0 is a diffeomorphism, we can rewrite above condition as < |(iet(^_)| < 00, where [A^l] = 
[dm^_/dx^\ , since A- = B- ■ d^Q. On the other hand, if Aq, is a graph over M ~ T^M, this can be 
further rewritten as < |det(^° )| < 00 , where [(A° )*^'] = [dm^_/dm^\ , since A^_ = A^ ■ A^^ and 
|(ief(^Q ^)| < 00. But since m_ = P-{t) = Expmi—v), we get (5.2), provided A^ satisfies (4.2) and r 
is given by the central map F(m) = v, generated by / via the central equation. □ 

Example: On ]R^, from (3.3) and (4.4), (5.2) becomes < |1 + (iet[5^/]| < 00, but since (5.1) must 
already be satisfied, we single out the new condition as det[d'^ f] 7^—1 , Vx G f7 C IP? . However, 
already on S"^ or H^, the explicit form of (5.2) in local coordinates, for generic /, becomes rather long 
and it is much simpler to check it directly, for each specific /, using the specific expressions for m_(m) 
obtained from (4.9) or (4.11). 

Conditions similar to (5.1) and (5.2) apply to nonstandard functions, as well. It is important to 

emphasize, once again, that any function satisfying the conditions of Lemma 5.2 does not necessarily 
generate a canonical transformation on M, since these conditions do not comprise a sufficient set, 
generically. A more complete analysis is needed for such characterization, which lies outside the scope 
of this work. Ultimately, though, we can check explicitly the consistency of the implicit map m+(m_) 
obtained from any function / on f7 C M, via the central equation and the symmetric exponential 
map. These functions shall be singled out : 

Definition 5.1 A function f = ^fa G C^(M), k > 2, which is (locally, on U C M) the central 
generating function of a canonical transformation on M is henceforth referred to simply as a (local) 
central action onM. The set of all such functions is denoted Ac{U). Specifically, A^{U) for standard 
actions. More generally, a function which locally generates a canonical relation on M , via the central 
equation, shall be called a (local) central relation on M, whose set is denoted TZelc{U) . 

Of course, AdUi) C AdU'i) if Ui D U2. Thus, VC/ C M, Ac{U) D Ac{M) the space of 
central actions on M, clearly non-empty since every function = —th, h G C|^(M), is a central 
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action, for t sufficiently small. On the other hand, since central relations must satisfy (5.1), but not 
necessarily (5.2), Ac{U) C TZelc{U). More generally, one might still wish to be free from constraints 
imposed by central caustics. For functions on M, this is generically impossible. On the other hand, 
inspired by some works in symplectic topology [|6|| [fzl] ||39[| , we could allow the generating functions to 
be defined on M x M'^, instead, using the extra variables to analyse the behaviour at those critical 
points. But the study of their general definitions and properties is not to be found here. In an 
independent context, a very important and particular case of such "extended" functions shall be seen 
in § 8,9, where d = 1. Back to functions on M, our main interest, in what follows, is concerned with 
their compositions. For central actions we need only worry about the presence of caustics, but for 
canonical relations we must also worry about the possibility of their compositions being well defined 



(clean products) |10|| |14||42||43|. Thus, for simplicity, we shall not pursue on this broader context here, 
focusing instead on the compositions of central actions properly, from now on. But we emphasize that, 
whenever well defined, the following rules of composition apply for central relations as well. 

Illustrations: The simplest examples of central actions are those which generate uniform translations 
on flat space. Let a = T^: IE? — > iR^, by x i— > which corresponds to the lagrangian plane v = ^/2, 

a constant, in TIB? . From (4.5) one has ^ = — J • (dfa/dx), integrated as fa{x) = (J • ^ ) ■ x , modulo 
constants, rewritten as a skew-product fa{x) = ^ A x . In particular (when = 0) the null action (or 
any constant) generates the identity, a fact valid for every M. The next simplest examples, still on 
flat spaces, are the homogeneous quadratic functions fa{x) = B x, where I? is a symmetric matrix, 
or equivalently fa{x) = I3p^ + f3'q^ + 27pg, where /3, /3', 7 € IR. When 7 = 0,^ = /?' = - tan(A/2), /„ 
generates a rotation by an angle A through the origin. Notice that /„, as well as (5.1), diverges 
everywhere when A = zbvr. This is a reflection at the origin, TZq, and the corresponding lagrangian 
submanifold in TIR? is At^^ = TqIR? which is not a graph over T^IR^ anywhere. When /? = /3' = 0, 
7 = — tanh(A/2), fa generates a pure hyperbolic transformation with stable/unstable submanifold 
coinciding with the p/q axis. Now (5.1) is always satisfied, but (5.2) fails asymptotically as |A| — > 00. 
It is not hard to see that every quadratic central action generates an element of the homogeneous 
symplectic group on IE? and, conversely, every such element which can be centrally generated, is 
done by a quadratic central action (Cayley transform). This is not a property of the more familiar 
generating funtions of mechanics on IE?, e.g. f{q-,q+), for which generic translations on M are also 
generated by quadratic functions. Adding the former two examples, fa{x) = ^f\x + x^Bx generate 
elements of the inhomogeneous symplectic group. On S"^, the simplest examples are the central actions 
for rotations. Here, standard case only. Let a = 0(^,27) be an element of the group of rotations 
50(3) acting on S"^ , whose pole (fixed point) is p and whose angle of rotation is 27. Taking polar 
coordinates for m = {9,f) and p = {x,^), and 7 G [— 7r/2, 7r/2], the central action is written as 

/„(m) = -25-1 {5^ [C\Cg + S^SgC^^^,)] } . (5.3) 

Although it looks complicated in local coordinates, f = ■^fa has the simple geometrical interpretation 
shown in Fig.5.1(a). If we notice that C^^Cg + S^SgC(^^^_^^ = Cy, y = distance{m,p), the map 
F„ : 52 ^ {TS\ is given by (4.7) as 

V = c-i lc,/^i-{s,Cyr} ^ C-^C./Cf} 1 

P = T-i {[S^CgC(^^_,) - C^Sg]/S^S(^^_,)} + (1 - Signiiip - e))7r/2 J 
and the implicit transformation m-|-(m_) generated by fa is given by (4.9) as 

0± = { [c^Cg ± s^SgS^s^^_,)] /^i-{S-yCyY] 1 

^± = ^± {5-, [C^Sg - 5;,CeC(^_,)] /A±} + (1 - Sign{\±))T: /2 , J 

where \± = C^Sg ib S^CgS^S(^^_^-^ . Notice that when m ^ p or p, v ^ and m is a fixed point of the 
transformation (5.5), as expected, V7 G [0, vr/2). However, when 7 = 7r/2 the transformation (5.5) 
is not well defined. Again, in this case q;(j?, 27) = a(p, vr) = TZp : S"^ — > 5^, and the corresponding 



(5.4) 



(5.5) 
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lagrangian submanifold in (TS'^)o is At^^ = (TpS'^)o, which is not a graph over T^S"^ anywhere. We 
can also see this singularity using (5.1), when 7 = 7r/2. To see that m ^ p\s not a caustic singularity, 
for 7 7^ '7r/2, one can check directly in (5.1). Similarly, expanding / around m in linear coordinates 
{PiQ)-, gives faifn) = constant + {T^)y'^ + o{y'^) ~ {T^)y'^, with y"^ = + q^, which is the central 
action for rotations on the plane. 




Fig. 5.1 - Full lines represent geodesies 



On H^, we now consider the central actions for some elements a G 50(2,1). A subclass of such 
elements is characterized by those a = a(p, 27), which are real rotations on around the real pole 
p G iJ^, through an angle 27 , 7 € [—1^/2 , '/r/2]. The corresponding central actions and canonical 
transformations are analogous to the spherical ones. But, in opposition to the real hyperbolic rotations, 
stand the ideal hyperbolic rotations. These can be characterized in two ways. We may consider 



a = 0(^,27) , 7 G [— 7r/2 , vr/2], p an ideal fixed point (a point in IR satisfying x + y 



while a real point in C ]R? satisfies 



+ y2) 



1, 



1). Alternativelly, we take a = d{i',e,z) 



where £ [0,oo), z £ (—00,00), e G [0,27r]; (z^, e) being the coordinates of a real axis in H , i.e. a 
geodesic in s.t. v is its distance to the origin and e is the angle this minimal geodesic arc makes 
at the origin, as shown in Fig. 5. 1(b). We interpret \z\ < 00 as a geodesic segment along this axis, its 
sign determining a particular orientation for it The corresponding central action is given by 



(5.6) 



where CySpC(^^_^) - SyCp = 5^ 



distance{m, axis). From (5.6) fa exists only when 



1^.5,1 <l^O<x<S-' {l/\S,\) , (5.7) 

but in this neighborhood of the axis, it is a well defined real function whose simple geometric inter- 
pretation is shown in Fig. 5. 1(b). From (4.10), the map Fa : —>■ TH^ is given by 



T-i {[C,Cpq^_,^ - Sjp]/C,S(^^_,^] + (1 + Signiif - e))7r/2 



(5.8) 
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and the /^-generated intrisic transformation m-|-(m_) is given by (4.11) as 



P± = c-^ 

ip± = (p±T 



A±} + (l + 5i5<±))vr/2 , 



(5.9) 



where ^± = CzSp ± SzCpCuS{^^^^) . The transformation (5.9) describes finite motions along lines 
equidistant to the axis (z^, e), as shown in Fig. 5. 1(b). We see that (5.9) has no real fixed point and is 
also free of central caustics, i.e. it is well defined everywhere on the same neighborhood of the axis 
{y, e) on which is well defined. In other words, the lagrangian submanifold C TH^ is a graph 
over this neighborhood; one can check explicitly that (5.1) does not diverge anywhere. Notice that, 
as X — > 0, m, m_|_ and m_ all lie on the axis, with z = distance (m,rn-4_). However, as \z\ — > oo, (5.2) 
fails asymptotically, in agreement with (5.7), which implies x = 0, when \z\ = oo. 



6. Composition of Central Actions 

Having explored the symplectic structure on central groupoids, we now add their algebraic structure 
in order to answer the following question: Let /c^, /^a £ C^(M) locally be the central actions for two 
canonical transformations, respectively, ai,a2 : M M, m i— > aj(m). Since a = 02(01) is another 
canonical transformations on M, what is, locally, its central action? In other words: How do central 
actions compose? We start by considering : 

Definition 6.1 Let Q he a (local) symplectic potential for uo on Uc_ <Z M . Then, on DU(^ C DM , 
6C := -P|(C) ~ P*iC) ^s a local symplectic potential for 6uj and := <I>q((^C) 'i'S locally a symplectic 
potential for f^o on C (TM)o, $o(^c) -^^C' ^'^ called an additive potential . 

The reason for this name is obvious once we realize that 

Z^Ti T2) ^ Z^n) + Z^(r2) . (6.1) 

Here, as in Def. 3.2, this is a shorthand notation for the fact that the local 1-form © 
vanishes on the graph of the groupoid composition. Such additive potentials do not coincide with the 
central potentials, i.e. does not satisfy property (4.1), even locally, and we cannot use it to define 
central generating functions directly. On the other hand, the central potentials are not additive, as in 
(6.1), but in order to compose central actions, we now realize that, on , 

Z^ = Zo + dQc , Qc^Q^^eC^niWO ■ (6.2) 

Therefore, if Tj G Aq,-, with Po^Ti) = mi £ U, is locally generated by the standard central action 
fai G C|j(M), then on {TU)o H W^: Zt^in) = dfa^nii) + dQt^in), combining (4.3) and (6.2). Hence, 
if Aq, satisfies (4.2) on (TU)o, for a = 02(01), and r = ri r2 € Aq, with Po{t) = m £ U, is locally 
generated by the standard cental action fa G C|j(M), then, on VF^ n {TU)q , from (6.1), 

dfa{m) dfa^mi) + dfa2{m2) + dxc{Ti,T2) , (6.3) 

where we define the standard phase function Xc = ^ (^f^)o ~^ locally by 

Xc(Ti,r2) = Qc(ti) + QdT2) - Qc(n T2) . (6.4) 

Now, we apply the following crucial result: 

Proposition 6.1 The function X( defined above (6.4) is well defined on the whole {TM)q and inde- 
pends on choices of local symplectic potentials on M. Actually it coincides, modulo constants, with 
the symplectic area of a standard geodesic triangle on M , i.e. which can be defined by elements in 
(TM)q , determined by its midpoints m,mi and m2. This area shall be denoted by Ao(?7i, mi, 771,2). 

Proof: Here we rely on a mathematical construction [Q] Q which consists of building an 5'^-principal 
fiber bundle over a symplectic manifold (M, w) with connection o whose curvature is w/A. We denote 
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this as {SM,a) (M,a;/A), da = tt*{oj/\). Such bundle is well defined only if f^x/B"^ ~ 

p G where B is any oriented, closed two-surface without boundary on M. Here, A is an auxiliary 
constant which can be set to zero in the end. To extend this construction to the double phase space 
[51], we identify SDM = {SM_ x SM+)/S^, {SM±,^a) being the same principal bundles over 
M, but for opposite connections, and the quotient is taken with respect to the diagonal action of 
C on SM_ X SM^. Actually, we have chosen the connection whose curvature is —5u!/\, so we 
should perhaps denote this bundle by SDM, but to simplify the notation we keep to SDM. Thus, 

S^ {SDM, [-5a]) H {DM,-5uj/\), d[-5a] = [tt]* {-5uj / X), whose elements are denoted by [a, a'], 
where (c, a') G SM- x 5M_|_. Then, choosing identity elements of the form [a,a\, we extend [51] the 
groupoid composition from DM to SDM by [cr, cr'j [cr', a"\ = [a, a"\. Now, we pull back SDM using 
the (restricted) symmetric exponential map to get a (trivial) bundle over {TM)q, denoted {STM)q = 

%{SDM) ""^^ (TM)o X SDM. Next, consider sections {TM)o {STM)o which are obtained via 
the parallel transport along the Po-fibers of identity elements in {STM)o, the latter being fixed by 
the above choice on SDM. These polarized sections £q : {TM)q (5TM)o are such that po{£q{t)) = 
[a', a"], where (a', a") are the endpoints of a horizontal lift in SM of the short geodesic from m' to 
m" , centered on m = Poir), for ^o{t) = {m',m"). Hence, if {ti,T2) € (TM)q, then eo(''"i) & £o{t2) = 
£o{tiQt2)Kq{ti,T2), where Kq : {TM)q S^ is the holonomy in SM over the triangle of composition 
on M. In other words, its symplectic area Ao{ti,T2), where the standard triangular area function A 
is a well defined continuous function on (TM)q. This means that Kq{ti,T2) = Exp{^^^Aq{ti,T2)} 
is well defined on the whole (TM)q , as well as being completely independent of local expressions for 
the connection [—6a\ and its pull-back a^ . But the pull-back connection ao can locally be written as 
ao — dO — jZ(^, where 9 is the fiber coordinate, for a local choice of symplectic potential u> 2± d(. Over 
the Po-fibers, these can locally be rewritten as ao\m — {dO — j{Z(^ — ^o))|m — {dO — jdQc^)\m, thus, 
provided QcIt^m = for a local representation of identity elements as (m,0;0), the trivializing 

sections Eq can locally be written as eo(''") — (t'i ^^Pi^^^QciT)}) and therefore the holonomy phase 
Aq(ti,T2) is locally identified as X^(ri,r2) = <5^(ti) -|- Q<^{t2) — Q^(ti T2). On the other hand, 
denoting a point in {TM)q by (mi, Vi; m2, V2), with T1QT2 = ts = (ma, U3), we can almost everywhere 
eliminate the three vectors by Expmii—vi) = Exp„i-^{—V3) = a, Expmx{vi) = Expm2{—V2) = b, 

Expm2{v2) = Expmsivs) = C, and ^o{mi,Vi) = {a,nm2'Tlmz{a) = h), ^>o(m2,V2) = {h,TZrn:,'Tlnii{b) = 

c), <&o(?Ti3, U3) = (a, TZm2'T^mi (fl) = c), the exception being those midpoint triplets //123 = {mi,m2,mz) 
for which the fixed point set of TZmiT^mzT^ms has nontrivial dimension. But for the set of such singular 
midpoint triplets, we can assume based on the lower dimensional examples, that the codimension of 
the set of singular midpoint triplets is always greater than 1. In other words, the map ^'o : {TM)q 
(M X M X M)midpoints is invertible on its immage set, except for a subset of codimension greater 
than 1. Hence, if /ii23 is a regular midpoint triplet, Ao(jUi23) = Ao(mi,m2,m3) = ^0(^0 ^(/^123))- 
Otherwise, if fii23 is a singular midpoint triplet, consider any continuous family of regular midpoint 
triplets A*'i23(e)) e > s.t. IJ^'usi^) ~ Mi23- These families always exist and, for any e 7^ 0, ^'q^ 
is well defined and continuous, so that Ao(ju'^23(^)) is * continuous function of e. It follows that 
Ao(Aii23) = Ao(mi,m,2,m3) = linic^o{Ao (//123(e))}, for any family /ii23(e) /U123 • Therefore, the 
holonomy can be identified with Exp{^^^^Ao{m,mi,m2)}, for nii = Po(Tj), m, = Po(ti T2), where 
Ao(m, mi, m2) stands for the symplectic area, modulo 27rA, of a standard geodesic triangle with given 
midpoints, for which XcC'^'Ij'^s) is thus a local expression. □ 

Let us remark that in the simplest case of iR^", the midpoint triangular area A is a unique well 
defined function on the whole iR^" x iR^" x iR^". Generically, however, Aq is well defined only on a 
subset U C M X M X M, as exemplified below for the torus, the sphere and the hyperbolic plane. 

In every case, though, for any triplet of points for which Aq is well defined, modulo constants, we 
get from propositions 4.1 and 6.1, via (6.3), the main result: 
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Theorem 6.1 Let fai, fai locally be standard central actions for a\, 012 M ^ M , respectively. If fa 
is locally the standard central action for a = 0:2 (ai), then 

fa{m) = /ai A/a2(m) = Stat(^^^^m2){fai{'mi) + /a2(m2) + Ao(m,mi,m2)} (6.5) 

defining the standard composition of central actions on M, if Ao{m,mi,m2) is well defined, up to 
constants, i.e. if {m,m\,m2) stand as midpoints of a standard geodesic triangle. 

Seen as a product on A%U), fai A fa2 may not exist or if so, it may not be unique, for C (TM)o 
may not satisfy (4.2) over U C M, or else it may be composed of many sheets in {TU)q. 

Furthermore, the composition rule (6.5) is not unique when considering nonstandard central actions 
as well, if M has closed geodesies. In this case, repeating the steps that led to (6.3), for reducible 
triples {i,j;k), using (4.3') we get the generalized version : 

d/^(m)«C^(mi) + d/4(m2) + d[xc]f,(fi,f2) , (6.3') 
where [xclij is a local representation for a general phase function {TM)fj — M, given by 

[xSjin,f2) = Qiih) + Qlif2) - Ql{c%{n,f2)) , (6.4') 

with each being a function defined as dQ^ = Z'^ — Zi, on appropriate subsets C {TM)i. As in 
the standard case, we identify this function by : 

Proposition 6.1' For (i,j;k) reducible, and for appropriate choices of Zi,Zj, and Z^ , the function 
bcc]ij ^^^^ defined on the whole {TM)f-, independs on local choices of ( and coincides, modulo 

constants, with the symplectic area of a reducible geodesic triangles on M which can be defined by 
elements in (TM)'fj x {TM)k and determined by its midpoints (mi, 7712; m), in other words, geodesic 
triangles with given midpoints and sides which are (i,j;k)-fold geodesies. Such area shall be denoted by 
A^^(m; mi, 777.2) or generically by A(m,mi,m2). 

Proof: This is an elaboration on the proof of Prop. 6.1 . Here, we must argue that there is a consistent 
choice for the Zj's, such that each Q-'^\j;jj^ = 7-' : M — iR is given by exp{^^^^7''(m)} = r]j{m), where 

r)j{m) is the horizontal transport over M 2± TiM of the Lj(mo) loop holonomy, given a reference 
point mo. Explicitly, if exp | \/— l/3j} denotes the Lj{mo) holonomy in SM (its homology class, if 
irreducible), then 7-'(m) = jY.L'{m) ™odulo 27rA and the constant fij , where dYjL'-{m) = L'-{m) = 

Lj(mo, m) is the reducible loop in M given as (mo — >■ m) o Lj{m) o (m — > mo) o Lj^{mQ), with Lj{m) 
denoting the j-fold geodesic loop based at m. Now, each <I>j : {TM)j DM provides a pull back 
(trivial) bundle {STM)j ^ SDM. Taking local trivializing sections ej : T^M {STM)j of the form 
£j = /37^([o", o"^]), where a'- = a ■ r]'j{m), for TT{a) = 7r{a') = m, which can locally be represented as 

£j{f^) — (m, \vj\] exp{^^^^^^ {m)}), then, from a similiar analysis to the standard case, the trivializing 

sections £j : {TM)j {STM)j which are obtained from ij via parallel transport along the Pj-fibers 

are such that Pj{£j{f)) = [a', a"], where {a', a") are the endpoints of a horizontal lift in SM of the 
j-fold geodesic from m' to m", centered on m = P.j{f), for ^j{f) = {m',m"). Hence, if (Ti,f2) G 
(fM)2^., then Cf^.(e,(fi),e,-(f2)) = ^^(^(fi, fs)) • i^i;-(ri, ^2), where C^^ : iSfM)^ ^ (^fM)^ j^the 
"extension" of the mixed composition c^j, defined as C^j{zi,Z2) := p^^(pi(zi) pj{z2)), Zi G {STM)i, 
etc... and K^-{f\,f2) € S^ is the holonomy in SM over the (i, j; /c)-triangle of composition, which, if 
{i,j;k) is reducible, corresponds to the symplectic area of this triangle. Since any reducible {i,j;k)- 
geodesic triangle stands in a bijective correspondence with a standard one, it follows that, as in the 
standard case, the (i, j; A;)-triangular area is a well defined function Afj of the midpoints. Therefore, 

the holonomy can be identified as Exp{^^^A'^j{m;mi,m2)}, for mi = Pj(fi), m2 = P7(t2) and 
m = Pfc(c^(fi, f2)), where A^(m; mi, m2) stands for the symplectic area, modulo 27rA, of the reducible 
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geodesic triangle with given midpoints and sides which are {k; i, j)-fold geodesies. This also means 
that K^j{fi,f2) = Exp{^^^A^j{m;mi,m2)} is well defined on the whole {TM)^- and independs on 
local expressions for the connection [—5a\ and its pull-backs oij. But, from the form of such pull-back 
connections and for Q(\j:j^j = T"' ^^'^ ^^^tIm ^ ^ above, it follows that the trivializing sections 
£j can locally be written as Sj{f) ~ (f; exp{-^^^(5;^(f)}) and therefore [x(]ij is a local representation 
for A^^(m; mi, m2), modulo 27rA. □ 

The appropriate choices of central potentials referred to in the statement are explained in its proof. 
Loosely speaking, each Zj is chosen as the pull-back of a single 1-form Z on TM, for each immersion 

{TM)j TM . Using (6.3') we get: 

Theorem 6.1' Let /^^, locally he central actions for ai,a2 '■ M — > M, respectively. If {i,j; k) is 
reducible and is locally a central action for a = 02 (oi), then it is given by 

ft{m) ^ fl,^ A'^fi^ = 5taV„^,){/^^(mi) + /4(m2) + Af^.(m; mi, ^2)} (6.5') 

defining the (i, j; k)- composition of local central actions on M, provided (mi,m2;m) can stand as 
midpoints for a reducible {i,j; k)-geodesic triangle and each action is defined with respect to the choices 
of central potentials referred to and explained above . 

Clearly, (6.5) is a particular case of (6.5'), with Aqq = Aq. However, by restricting to the standard 
composition we have obtained a uniquely defined rule. In extending to generic central actions this rule 
ceased to be unique, even though each instance is well defined. Thus, for simplicity, when considering 
multiple compositions, in § 7, we keep to the standard case. Finally we remind that, whenever 
meaningful (clean product |4^J, T^-^'j applies for the composition of central relations on M. 

Examples: Standard cases, only. In the simplest spaces, we provide the expressions for A and analyse 
the specific nontrivial restrictions which apply, in each case. First : 

A (x, xi, X2) = 2 {x A xi + xi A X2 + 3^2 A x} (6-6) 

is the midpoint triangular area on defining, via (6.5), the composition of two central actions on 
M^. On T^, we add the restriction on the lengths of the triangles, in terms of midpoints : 

ki - 9I , |Pi -Pl , ki - 92I , |pi -P2I < vr/2 , i = l,2 . (6.7) 

On the 2-sphere, we have (see Appendix ) : 

A = 2Arg |o-\/l - 1)2 + (6.8) 

where a = a{m, mi, m2) G {—1,0, 1} has the same sign as any of the scalar products (m-mi), (m-m2), 
(mi • 7712), and D = D{m,mi,m2) = det[m,mi,m2]. Here, rh £ S'^ C M^. For small triangles, i.e. 
|A| < TT, (6.8) simplifies to A(m,mi,m2) = 25'""'^ (det [m, mi, 7712]). This explicit form for A must be 
placed in (6.5), provided the triplet (m,mi,m2) stand as midpoints of a spherical triangle with all 
sides shorter than tt. This nontrivial restriction : 

sign{{m ■ rhi)) = sign[{rhi ■ 7772)) = sign{{m2 ■ rh)) (6-9) 

is the constraint that must be imposed on (6.5) and (6.8) to define the standard composition of two 
central actions on the sphere. Constraints analogous to (6.9) or (6.7) apply to each mixed composition, 
for every space with closed geodesies. The noncompact hyperbolic plane, though similar to the sphere, 
presents one subtle distinction. Here we have, for m ~ tti G C IB? , 

A(m, mi, m2) = 25""*^ (det [777,, mi, 7772]) . (6.10) 

But now we notice that (6.10) only makes sense on the subset of x x determined by : 

I det [777, 7771, 777-2] I < 1 , (6-11) 
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The compostion of two central actions on is given by (6.5) and (6.10), subject to constraint 
(6.11). Actually, (6.11) is a necessary and sufficient condition any triplet of points on must satisfy 
in order to be the set of midpoints of a hyperbolic triangle. As in the case of the sphere, there is a 
nontrivial constraint in the composition of two central actions, but in the hyperbolic case, (6.11) comes 
naturally from the definition of the area function (6.10), being intrinsic to the hyperbolic geometry. 
Its interpretation is that the three midpoints can never be too far apart, as can be pictured by drawing 
a geodesic triangle on the Poincare disc. For the sphere, on the other hand, (6.9) is a consequence 
of the restricted groupoid (TS'^)o. In other words, there do exist spherical triangles whose midpoints 
do not satisfy (6.9), but they cannot be used in the definition of the standard composition of central 
actions (they can be used for mixed compositions), contrary to the hyperbolic case where, for every 
triangle, (6.11) holds. Also note that, when defining triangles by their triplets of vertices, instead of 
midpoints, no non-trivial restriction exists (except for sets of measure zero), which is a consequence 
of the fact that the map {vertices} — > {midpoints} is generally less trivial than one would guess by 
looking only at the euclidean plane. See also Appendix and []3l|p3|[38|[47|. 



7. Multiple Compositions 

So far we have seen how two central actions may compose into a new one via (6.5'). Now, we 
want to generalize it for the composition of an arbitrary number of standard central actions, whenever 
possible. This should be done with some care, so we first get : 

Lemma 7.1 Let fa,, locally be standard central actions for Aq- C {TU)q. If fai A and A fa^ 
are unique standard central actions on U, as well as {fai A /aj) A fa^ and fa^ A (/qj A fas), then 
{fai ^fa-i)^ fa-i = fai A {fa^ A fas) is given by 

fai A fa2 A fas (m) = Stat(^ 

mi ,m2 ,ms) 

\^faA^i) + P4im,mi,m2,m3)\ , (7.1) 

Pi being the smallest symplectic area of any ( oriented) quadrilateral decomposable in triangles defined 
by elements in (TM)q only and determined by the midpoints (m, rn-i, m2, 772-3), up to constants. 

Proof: Suppose fai A/qj, given by (6.5), is a unique standard central action on [/ C M, just as fa^- 
Applying (6.5) again: {fai A fa^) ^ fasim) = Stat^^j^, ^^^){fai A fa^im') + fas{m^)+ A{m,m',m3)} = 
Stat(^rn' ,ms){Stat(^^^^^^) {fai{mi) + fa2{m2) + A(m', mi , 7712)} + /as ("is) + A(m', ma , m)} . If a unique 
solution exists, this rewrites as: 5tat(„j_„2,m3){/ai("T'i) + /a2("^2) + /o3(^n'3) + 'S'iai(m'){A(m',mi,m2) 
+A(m', ma, m)}}. But, with {rrii} constrained by the overall stationary condition then, via the central 
equation (4.3), A(m',mi,m2) =: gi{m') and A(m',m3,m) =: g2{m') provide well defined maps Gi, 
G2 ■ m' iji,V2 £ (Tm'M)o, supposing that both gi and 52 , as well as fai A /q-j and {fai A /^j) A fa^ 
are well defined unique central actions, i.e. both partial and complete unique solutions to the stationary 
conditions exist. From the form of the central equation and the involutive character of the central 
potential: Zq{t) = —ZQ{f), the m' stationary condition implies Gi(m') = —G2{m'). Using the 
symmetric exponential map, we see that the two triangles composed form a single quadrilateral, i.e. 
their sides centered at m' are precisely opposite to each other. In other words, we have that 

Stat(^m') { A(m', mi, m2) + A(m', ms, m)| = P4(mi, m2, ms, m) (7.2) 

where P4^{mi,m2,rn3,m) is the symplectic area of a standard quadrilateral whose midpoints are 
(mi, m2, ms, m). Generically, the four midpoints do not determine the quadrilateral uniquely. How- 
ever, if all of its triangular decompositions yield triangles which are defined by elements in (TM)q 
only, then by Proposition 6.1 and (7.2), P4 is the symplectic area of any such quadrilateral, up to 
constants : any degeneracy in its specific geometry, which is a continuous function of m', does not 
alter the symplectic area, and any other "(TM)o-quadrilateral" , in the above sense, will have the same 
symplectic area up to constants. □ 
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The existence and uniqueness requirements are quite difficult to assure beforehand, in general. 
Of course, if any /„. or intermediary composition is not a central action, the triple composition is 
void. On the other hand, if an intermediary composition is not unique, say fai A = {51,52}) 
then we could proceed to {gi A fas, 92 A fas} in just the same manner, but the final composition 
{fai A /eg) A fas is not guaranteed to be associative, in principle. However, if all A^,. , Aq.q,^.^... are 
composed of sheets whose one of them can be consistently singled out, for some particular reason, as 
well as their corresponding central actions, then we can apply Lemma 7.1 exclusively to this particular 
set. Carefully reiterating all steps to (7.1) gives : 

Corollary 7.1 Let fai locally be standard central actions for Aai C {TU)q. If all intermediary, ordered 
compositions fai A fa,+^, fai A {fa,+i A /o^+a), {{fai A A /a.+a) A fa,+s' ^^c... are unique standard 

central actions on U , then any n-string of ordered compositions which is a unique standard central 
action equals any other such n-string and is given by: 



/ai A A ••• A /Q„(m) = 5tat({^,j) |^/Q,(mj) + P„+i(m, {mj) j (7.3) 

Pn+i{m, {mi}) being the smallest symplectic area of any (oriented) {n + l)-polygon which can be 
triangulated by elements defined in (TM)q only and determined by the midpoints {m, {rtii}). 

Again, if uniqueness fails, but a unique set of central actions can be consistently singled out, then 
we can use corollary 7.1 for these particular standard central actions exclusively (see §8). We should 
note that, with greater care, the above rules can be generalized to multiple mixed compositions of 
central actions and, even more generally, central relations whenever meaningful. 

Corollary 7.1 generalizes a previous result on iR^" |27]. We remark that some of the above discussion 



regarding composition of midpoint triangles, in more general symmetric symplectic spaces, has been 



approached independently from the point of view of associativity for star products [31| 



Examples: For the composition of three actions, the euclidean plane presents an interesting fea- 
ture. From (6.6) and (7.2), we get P4/2 = Stat(^/) {xi A X2 + X3 A x + x' A {xi — X2 + X3 — x)}, which 
implies: xi — X2 + — x = . This means that (x, {xi}) are the vertices of a parallelogram 
with diagonals X3 — xi, x — X2- But this is true for any quadrilateral in M^, i.e. their mid- 
points are vertices of a parallelogram. Conversely, given any parallelogram on M^, there exists 
a continuous family of circumscribed quadrilaterals whose midpoints are the vertices of the given 
parallelogram. Such a family can be parametrized by one of the vertices of each circumscribed 
quadrilateral or, equivalently, by the midpoint of one of its diagonals, as x' above. In accordance 
with (7.2), the symplectic area independs on x' , being uniquely given as twice the area of the in- 
scribed parallelogram: P4{x, xi, X2, x^) = 2{xi A X2 + X3 A x = x A ifi -|- X2 A X3), which, inserted 
in (7.1), defines the composition of three central actions on M^, with {x,{xi}) subject to the paral- 
lelogram relation and under the necessary existence and uniqueness conditions presented in Lemma 
7.1. On the torus, the same analysis and results apply, but now subject to the extra (standard 
groupoid) constraint: \q — qi\,\qi — qj\,\p — Pi\,\Pi ~ Pj\ < ^/^ > hj = 1,2,3. In the spherical 
case, on the other hand, the above quadrilateral ambiguity is an exception, when considering only 
those quadrilaterals defined by composing elements in (TS'^)o. Then, as with spherical triangles, 
a near-bijection between the sets of vertices and midpoints allows for a unique definition of the 
quadrilateral geometry either way, i.e. each quadrilateral is uniquely determined by its vertices or 
its midpoints (with a few exceptions) and the only restrictions on the latter derive from the re- 
stricted groupoid (rS'^)o. Denoting Cij = rhi ■ rhj = cosine{distance{mi,mj)), we can write the 
midpoint area of convex quadrilaterals which are decomposable in small triangles, i.e. |A| < vr, as 
P4(mi,m2,m3,m4) = 2cri234 • C'~^{Ci2C34 + (723(^41 - C13C24}, where cri234 = ±1 is the orienta- 
tion of (mi, m2, m3, 7714) and the standard groupoid restrictions on the midpoints of these simplest 
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quadrilaterals now become: -D123 , -D234 , -D341 , -D412 > , where Dijk = det[rhi,mj,fnk]-, which, to- 
gether with (7.1), define this simplest composition of three central actions on S"^, under the existence 
and uniqueness assumptions of Lemma 7.1. On , the midpoint area function for a convex quadri- 
lateral is given similarly by Piimi, 1712,^13, 1714) = 2(Ti234 • C~^{C'i2C'34 + C23C41 — C'i3C'24} , where 
Cij = cosh{distance{mi,mj)). Again, the restrictions are intrinsic to the hyperbolic geometry and 
can be obtained directly from the area function, i.e. the convex set (777-1,7712,^3,7714) must satisfy 
|C'i2C'34 + ^^23^^41 — C'i3C'24| < 1 and each quadruplet of midpoints satisfying this constraint deter- 
mines a unique convex hyperbolic quadrilateral and vice- versa (with a few exceptions, see below). 
Using Lemma 7.1, we obtain this simplest composition of three central actions on H^. See Appendix 
for a more detailed analysis on these quadrilateral geometries. 

Besides providing explicit equations for compositions of three central actions, the previous discus- 
sions further illustrate some kinds of constraints which the midpoints, or centers (the arguments of 
the composing actions) are subject to. In this respect, the euclidean plane presents the feature that, 
when the number of composing actions is even, no restrictions apply, but when the number is odd, 
there is a degeneracy in the determination of the 2£-polygon from its midpoints, corresponding to a 
linear functional restriction g{x, xi, • • • X2i-i) = on the arguments of the composing central actions. 
For the torus, one must further add the groupoid restrictions. On the other hand, in the nonflat 
cases studied, such a degeneracy is an exception and its corresponding extra constraint is not present. 
When considering only those polygons which can be obtained by iterated (standard) central groupoid 
compositions, there is a near-bijection between the sets of midpoints and vertices. On H^, intrin- 
sic restrictions on the midpoints exist, though, which show explicitly in the midpoint area function, 
whatever the number of composing actions. On S"^, only the (standard) central groupoid restrictions 
apply, in every case. 

Another way to view degeneracies and extra constraints is the following: Let {m,{mi}n) be a 
candidate for the set of midpoints for an (77 + l)-polygon on M. Also, for each 777j, let IZrm ■ M ^ M 
be the corresponding involution whose fixed point is 777,. Then, the existence of a circumscribed 
(77 + l)-polygon to the midpoints (777, {777j}„) is equivalent to the existence of a fixed point for the 
symplectomorphism Vn+i = T^mi • T^m2 ^m„ • T^m '■ M — > M. Now, on M^, when 77 is odd Vn+i 
is a translation, see (2.1). Fixed points exist only when this translation is the identity, in which case 
every point is fixed. For the sphere, on the other hand, Vn+i S SO{3), V77 G IN, and there is always 
a fixed point (actually two). However, for H^, Vn+i S 5*0(2,1) and there may or may not exist a 
fixed point on H^, \/n > 2, but when there exists, it is unique. The exception, in both cases, is when 
Vn+i = 1 and for these sets of points the corresponding circumscribed polygon is not uniquely defined 
by the midpoints, but so is its area. 

We shall not present here an explicit characterization of multiple compositions, for generic 77, in 
every example. The reader is referred to for the euclidean case. Instead, in the next paragraph 
we study a particular limit for (7.3), as 77 ^ 00. 

8. The Central Variational Principle 

We now focus on the relationship between finite and infinitesimal canonical transformations in the 
central description, i.e., on the relation between finite and infinitesimal central actions. To this end, 
consider (a)r = {q^^^^It; a continuous sequence of canonical transformations on M for which the 
following properties hold: Vx, y,z,t ^ [0, T] s.t. x + y + z = t , a^^^ : M ^ M satisfies {a'y'^)*uj = u>, 
at = a[^^ = ai^^^"* (oy^'' • a'^^) = a'z^^'^ ■ <^xXy = [Gt'z^^'^ • alf"*) • ot'x' = ot'y+z ' '^^'^ and also at+e = 

a^P ■ at = at + o(e), as e ^ 0, with ai*"* id : M ^ M, where a'^^ is the infinitesimal canonical 
transformation defined by the hamiltonian function h{t), via Hamilton's equation or, equivalently, by 
the infinitesimal central action f^(t) = —eh{t), via the central equation (4.3). We assume that h{t) 
is a continuous function of t, but h{t) 7^ h(t'), in general. In other words, /i is a nonautonomous 
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hamiltonian, i.e. a continuous function on M x [0,T]. Accordingly, we denote /i(0) = h and, if 
h{t) = h, yt ^ [Oj^li we say that h is autonomous. In this particular case, h{t) = h, we have that, 
for t £ [— T, r], {at)~^ = OL_i. But generally, i.e. 7^ /i(i')> such simplest inversibility relation only 
applies for very short intervals of time, i.e. (ai^^'^ = q^^ only in the limit e — > 0. In other words, the 
flow of {at}T is locally hamiltonian (autonomous), but not globally. Now, we will seek local central 
actions for at given in terms of h{t). To achieve this goal we shall use the results from the previous 
paragraph, but, in doing so, we should certify that fa^ exists and is unique, \/t' G [0,t), in principle. 
Actually, both conditions can be relaxed in this particular case, as is shown below. Then, we obtain: 

Theorem 8.1 Let {a)T be a continuous sequence of canonical transformations on M , as above, where 

as e ^ 0, af^ is the infinitesimal canonical transformation generated by the hamiltonian h{t) = h{m, t) 

continuous in t. Then, wherever the central action for at = af'^ exists, fat{i^) = = ^hi^^t); 

it satisfies the Central Variational Principle: 

'^h{m,t) = Statf^^) i^- j h{m'{t'),t')dt' + ^uj^ {m,t) , (8.1) 

for a family of continuous paths 1/ : [0, t] — > M geodesically centered on m, where, by definition 
{i^uLo}{m,t) = [i^,^ujY(m) is the symplectic area between the curve v and the geodesic from v{t) to z/(0) 
centered on m. This area function is well defined, up to constants, provided the geodesic is such that 
the full closed circuit reducible. Furthermore, the stationary paths v solving (8.1) coincide with the 
classical trajectories on phase space describing the continuous evolution from z^(0) to u{t) . 




Fig. 8.1 - Central symplectic area of a curve 



Proof: We start by considering those central actions which are continuously deformed from the hamil- 
tonian, i.e. fat — *■ —sh, as t ~ e — > 0, continuously on [/ C M. Thus, let {Aq^} be a continuous 
set of lagrangian submanifolds in {(TM)q,Qq), evolving from the zero section T^M = A^g ~ M, s.t. 
Vr G Aq,(, ^{t) = (m_,m+) satisfies m+ = Q!f(m_). Suppose that over U C M there are £ branches 
of A„ = Aq,^, where a„ = at for t = tn & (0,T], generated by i central actions fi„,j = I-- - ,i. 
Generically, U is an open neighborhood bounded by the sets of central caustics j = 1, • • • £ — 1. 

Further, let t^ £ (0, t„) be the time when a central catastrophe first appear for {Aq^}, i.e. A^t satisfies 
(4.2) everywhere on M, for t < t^ only. Denote {m^} the set of central caustics for A^, then we can 
find neighborhoods in M divided by {rrik}. This is not so clear when codim({mfc}) > 1, but following 
the evolution of these caustics we can find appropriate subdivisions in M. Let's denote by C/^ one of 
such neighborhoods, containing U . Then A^^{U) D A^^{U^). Let /^^ G •^c(^) be the central action for 
a^ on U d C M . Now, of all the branches of A„, over U , only one is a continuous deformation 
from A/j, over U , therefore, of all the central actions for a„ on U , only one is a continuous deformation 
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from /^^ on U. And since /^^ is continuously deformed from —sh, even though there are more than 
one central actions for a„ on U, there is a unique central action for on U, denoted /^^, which 
is continuously deformed from /^^ = —eh, t ~ e ^ 0. This result clearly independs on the number 
of subsets of A„ over U, or equivalently, on the number of central actions for a„ on U. Also, the 
other sets of central catastrophes for A„, {r„}-^, j = 2,---£ — l, are not connected to A^, so /^^ is 
actually a continuously deformed central action on the larger neighborhood z:) U, C C M 
being limited only by the set of central caustics {mj/}^, Vt' G (0,tn). This also means that, for such 
values of t' , /^^, exists, besides being continuously deformed from the hamiltonian, on U^^ C M. 

We are now ready to apply corollary 7.1 to = f^(r-i) A f^(r-2)' ' ' ^ f^a) ^ f^(o) i ^ times. That 

a)' a;, a;/ a) / 

tjr tjr tjr tfr 

is, the continuity of has substituted for uniqueness. From the involutive character of the cen- 
tral potential Zq, via the central equation we have that /(ae)-i = /a-s = ~fas, which implies that 
fa_^ = +sh + o(£^). Also for / ^t') = —eh{t') + o{e^). Then, letting r — oo, we can safely make the 

approximation /\,) = -^/i(t-) + o{{t/r)^) ^ -^/i(t-)> get on C M : 



tjr 



r-^oo 




Stat({^/}^) \2_^\ — ) h{mi,t'i) + Pr+i{m,{m'i}r 



(8.2) 



where m'^ = m'{t'j), t\ < t'^j^^ G [0, t]. Now we realize that (8.2) takes on the form (8.1) as we identify 
V = Imxr^ooiim' {t[)}r) '■ [0,t] ^ [// C M as a continuous curve satisfying Po{^Q^{i'{0),i'{t))) = m. 
In passing from (8.2) to (8.1), this area is integrated via a limit (r oo) of the area of (r + 1)- 
polygons whose midpoints are (m, {m'-},.), as r of the sides tend to length zero while the other tends 
to the geodesic from to m'l centered on m. In order to see that the stationary path is the classical 
trajectory with endpoints geodesically centered on m, we notice that each small side of the (r + 1)- 
polygon is a geodesic which, in the limit of very short time intervals, i.e. of very small sides, coincides 
with the local hamiltonian flow of h{m\,t\), which is centered on as discussed in § 4. Thus, 
taking the limit r — >■ oo of {m'(t'j)}r, we obtain a path that is everywhere tangent to the locally 
hamiltonian flow, in other words, that converges onto the classical trajectory. At first, (8.1)-(8.2) 
would apply only to those central actions that can be continuously deformed from the hamiltonian 
h, on U} C M. However, if fa^ is not of this type, it is always possible to decompose it in the form 
fat = /^rt'i ^f^(o) or f^(t'+t") if^(t') ^f^(o)) ,etc... where each "smaller" component is continuously 

deformed from the appropriate h{t'), on each appropriate neighborhood. Thus, they can be written 
as solutions to (8.1), with each z/' : [t',t — t'] — ^ M, etc... However, the laws for compositions of 
central actions (6.5'), plus the fact that we are composing central actions for the continuous sequence 
(q;)t, imply that the trajectories ly',^" compose into a single continuous trajectory u = v" o u' , for 
appropriate choices of fa,jf^(t')j with the corresponding areas and ^"W summing up to ^w. 

To see this, we notice that the stationary condition on m' in (6.5') implies, via central equation, 
that [(^w]* (m') and A(m',m",m) provide maps from m' into reciprocally inverse elements in T^^'M, 
which, via the symmetric exponential map, tells us that they compose into a single geometric figure. 
Repeating the analysis with m" , we have that (m'), [^u;]* (m") and A{m,m' ,m") compose 

into [^w]*(m) = {•^i,w}{m,t), for t = t' + t", using the fact that u = oi/', provided the composition 
is reducible, i.e. provided the triangle of composition A = A^j is a reducible circuit, which means that 
the geodesic from iy{t) to iy{0) centered on m is such that it closes the trajectory v into a reducible 
circuit. And so on... it follows that fat can also be written as a solution to (8.1), even when it is not 
continuously deformed from the hamiltonian, meaning that, for some {m',t') G M x [0,t), fa^ii,^') 
does not exist, i.e. m' is a central caustic singularity for Aq,^, , t' < t. So, even though generically /q^, 
does not exist everywhere on M , for t' < t, and f^t is not everywhere deformed from the hamiltonian 
function h continuously, wherever /^^ exists it can be written as a solution to (8.1), with v being a 
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continuous trajectory, v : [0, t] M, and m being the center of the (short or long) geodesic from v{t^ 
to 1^(0) closing the trajectory into a reducible circuit. □ 

Theorem 8.1 generalizes to nonflat symmetric symplectic spaces the previous result on euclidean 
space |2^. It is a real variational principle which is invariant at least under general transformations 
on M preserving the affine connection and the symplectic form. Besides, it does not require any 
local decomposition of the phase space M into lagrangian subsets. That is, it is fully adapted to the 
nontrivial geometry of M. Furthermore, in opposition to the complex conterparts, this real variational 
principle has only real classical trajectories as stationary solutions. The novel feature is that such 
trajectories are constrained on their geodesic centers instead of the more familiar (local) lagrangian 
coordinates of their endpoints. Therefore, in solving for the paths v which are stationary in (8.1), only 
the time t and the center m are held fixed. Finally, the central action ^h{m,t) itself provides, via 
the central equation, the finite transformation j/(0) — > v{t) and, given its explicit relationship to the 
infinitesimal generators h{m,t), plus the fact that it is a real function on M x [0, T], can be seen 
as a finite time extension of the hamiltonian function. 

9. Temporal Evolution of Central Actions 

We have just seen how the central variational principle provides, not only for the classical trajectories 
obtained by the stationary condition, but also for the central actions which generate finite canonical 
transformations and can thus be seen as finite time extensions of the hamiltonian functions. Now, we 
shall investigate the temporal evolution of such central actions. This can be done in two ways. First, 
we can examine the temporal evolution of ^h{m,t), for fixed m. The total variation of ^h{m,t) with 
respect to t depends on the direction of m. We have : 

e-Vt{^h{m,t))=e-\d^h{m,t)/dt + i^\d¥f,{m)^ , (9.1) 

denoting 5t = e . On the other hand, by (8.1), 

e-Vt{^i{m)) = -e ■ h{u{t),t) + e ■ \Th\d[.^w\\m)^ + o{e^) , (9.2) 

since ^^(m) is stationary in v and only terms in {Sv)"^ contribute, where o (((5i^)^) ~ o(e^). In equation 
(9.2) we are thus approximating the new path by the old one vt , i.e. we consider only infinitesimal 
variations in the endpoint u{t) along the same classical trajectory z/ : (0, t + e) ^ M. Via the central 
equation, we have = dfat : m i— > f^j, while dl^^wY : m ^ f^. But, by construction, = = 
$~^(z^(0), i/(t)) and therefore d^'^(m) = dl^^i.w]^ (m) , where now we identify ^'^(m) = ^\{m,t) as a 
generic central action. From this and (9.1)-(9.2), we identify: 

d'i'h{m,t)/dt = -h{m+,t) (9.3) 

But since m+ = Expmi+F^^{m)) , F^^{m) £ {TmM)i defined by d/^^ = via central equation, 

H [^l(m, t)] := h [Expm (m)) , t) (9.4) 

defines the functional M on Ac{U). Then, (9.3) can be rewritten as: 

d^h{m,t)/dt + H[^h{m,t)]={) , (9.5) 
which is the central version of the Hamilton-Jacobi equation . 

Examples: On (9.5) becomes d^h{x,t)/dt + h{x - \ J ■ {d'^h{x,t)/dx),t) = 0, see |2l||^. On 
5^ and H^, however, its explicit generic form in local coordinates is quite complicated and it is rather 
simpler to use (4.9) or (4.11) to write lH[^j^(m)], for each specific h. As simplest example, consider 
on S"^ the hamiltonian function h = —Cq, generator of infinitesimal rotations around the south/north 
axis, or poles. Using the convention \I'^ = f^^ = 2/j, by (4.9) we rewrite the standard version of (9.5) 

as df{e,t)/dt - i CeJl- {df{e,t)/def = , where we denoted ft{0,ip) = f{0,t), exploiting the 
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(/3-invariance of the action; and remembering that, by rescaling ^'^ = /q.^, we must also rescale h, and 
hence IH by the same factor. Check exphcitly that the central action given by (5.3), with x = ^-nd 
t = 27, satisfies the previous equation. Similarly for generic examples. 

Finally, we notice that we could rewrite equation (9.1) as: Vf (^'^(m)) := (9^^(m)/(9t + -UgJ(i^J^(m), 
for m = Vg defined by Vg\w = —dg, g G C^{M). Using (9.5), we get: 

Proposition 9.1 Let g,h £ C^{M x [0,r]) and ^{ £ Ac(U), U C M, where ^{ = ^tit) is related 
to h via the central variational principle (8.1), for a given time t £ [0, T]. The time derivative of 
"along g", i.e. in the direction of the local hamiltonian flow of g, denoted (^'^), is given by 

VU^hit)) = {^h{t),g{t)} - Mi^^M , (9.6) 

where { ,} denotes the Poisson bracket and the functional M is defined by (9.4) via the central equation. 

The r.h.s. of (9.6) involves partial derivatives in M only. Particularly interesting is the case where ^'^ 
is known to be invariant in a given direction Vg , 

Vfi^hit)) = 0, giving {■i'h{t),gm = M[^h{t)] (9.7) 

as a direct relationship between the functional JH and the Poisson bracket with g. Conversely, any 
function g satisfying the second part of (9.7) defines curves Tg : [0,T] ^ U C M , along which ^'^ is 
constant. The other particularly interesting case is when {'ifh{t),g{t)} = 0, for which Vf (^'?j(m, t)) = 
9^'/i(m, t)/5t, computable by the Hamilton-Jacobi equation. 

10. Conclusion 

We presented a general formalism for describing hamiltonian systems defined on symmetric sym- 
plectic spaces, where the local generating functions are real functions on phase space. We saw how 
the central actions (relations) are defined, generating finite canonical transformations (relations) via 
the central equation and the symmetric exponential map, and how they compose via a neat formula 
involving the midpoint triangular (polygonal) area. Then we saw how the "extended hamiltonians" 
satisfy a geometrically simple real variational principle, which determines the classical trajectories, 
and obey a Hamilton-Jacobi equation, mixed with Poisson brackets. 

The authors' main motivation into this central formalism lies in its application to problems in quan- 
tization and semiclassical analysis, more specifically in connection to "Weyl quantization" and "star 
products", which attempt to understand the classical-quantum relationship within the phase space 
formalism. In this respect, and specially for oscillatory phenomena, it turns into a definite advantage 
the definition of real phase space generating functions, which can be connected to hamiltonians in 
such geometrical fashion, with their neat triangular law of composition (see fsll [p7| pl| [p2| [^] [p7| ) . 

A possible extra application refers to implementing new symplectic integrators (see for a review). 
Here, the polygonal law of composition (7.3) could be applied in the discretization process, making use 
of the fact that each local action for a finite small interval of time is a small, in principle controlled, 
deformation of the hamiltonian. Also, for autonomous hamiltonians the formalism is symmetric with 
respect to trajectories in both temporal directions. 

On its own, however, such geometrically simple law of composition presents new routes of investi- 
gation on symplectic dynamical systems. For, as hamiltonian functions correspond to generators of 
infinitesimal canonical transformations, the actions correspond to elements in the Lie group. There- 
fore, their homogeneous presentation, which naturally extends to relations, sheds new light on the 
canonical formalism on symmetric symplectic spaces. 
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Appendix 



Triangular computations : In the spherical case, recall that if Aj are the angles of a geodesic triangle 
on 5^, whose opposite sides are £i < it, respectively, then the following trigonometric equalities hold: 
SiJSx^ = constant, Ci^ = Cefie^ + Se^Se^Cx^ , Cx^^ = Sx.Sx^Ci^ - Cxfixy Also, the area of a 
spherical triangle is its excess angle [36]. Now, for any spherical triangle with sides £i < vr, a simple 
computation shows that T^. = K/ (Q. —Ci^Ci^) where Aj is the angle opposite to £j and K isa. constant 
for this triangle, = 1 — — Cj^ — + 2(7^^ C^jC^g = Det^[ai, 02, 03], where are the vertices of 
the triangle. By correctly fixing the orientations, we can take the +sign in the square root. Then, let's 
denote li = 2yi and Xi = distance{mj , mk) , where mi is the midpoint of the £i side. We have that rhi = 
2n—i^j + "fc), from which, since C^^ = rhj ■ rhj,, we get that Cx, = 2^i^y, + ^"2,) = 2U~(^yk + C'^J, 
where zj = distance{mj,aj). From this we get: Cxi/Cy-^ = C^j/Cyj = Cx^/Cy.^ = F, a constant . 
This is a generalization of the plane trigonometric relation Xi/yi = 1. Now, in order to compute F, 
we substitute the previous equation in the trigonometric equalities, to obtain: F^ = C^^ + C^^ + 
Cx^ — 2CxiCx2Cx3 = 1 — -Det^[mi, m2, ma]. Using the previous equations for Tx,, F and F^, we finally 
get Ta, = F^/r^/(F2 - Cx.CxJCx,) and Ti^x,+x,+x,) = Ta = vVT^/iT'' - 1/2). Identifying 
F = C^, we have Ta = T2^, that is Ca/2 = ±\/C'xi + C'x^ + C'^^ - 20x^0 X2CX3 ^ F , 5'a/2 = 
Det[rhi,7fi2-,rh^- The sign choice for S'a/2 is fixed by the orientation. We still have to determine the 
sign of the square root in Ca/2- Obviously, if |A| < vr, we must choose the +sign. These triangular 
areas are continuously deformed from infinitesimal triangles, for which Cx^ > 0, Vi. Since we are 
considering only short triangles, i.e. yi < tt/2, Vi, from CxJCy. = F, we get F > 0. On the other 
hand, let |A| = 27r, i.e. consider mi,m2,m3 to be collinear, same for ai, 02, as, such that the 
"triangle" coincides with a geodesic meridian. Again, if yj < 7r/2, Vi, it is clear that in this case 
Cxi < 0, Vi, from CxJCy. = F, since Ca/2 = F < 0. And so on for triangular areas continuously 
deformed from this "big triangle". Finally, when |A| = vr, Ca/2 = and we have that, Vi, Cx^ = 0. In 
this case, Cy. is completely undetermined and so is the triangle, although the area of all such triangles 
is uniquely given by their common midpoints. It follows that the sign of the square root is the same 
as the sign of the Cx^, or in other words rhj ■ rhk, '^i,j, k, if all sides are short. Hence, we've got (6.8). 
As for the restrictions on the midpoints, we have already obtained these. If Cy- > 0, Vi, then either 
C^. > 0, Vi, or Ca;. = 0, Vi or else C^. < 0, Vi. In all cases, sign{CxJ = sign{Cxj), yi,j, which is 
condition (6.9). With some care, these analysis and results can be extended and modified for general 
spherical triangles. Again, this previous analysis can be adapted to H^, with care (refer to [^] for 
hyperbolic geometry and trigonometry). On we don't have problems of antipodals or sign choices, 
since |A| < vr, always. On the other hand, the analogous to (6.8), namely (6.10), is well defined only 
when (6.11) is satisfied. But once F = Ca/2 = C^^ + C^j + C^g — 2Cxj^Cx2Cx^ is well defined, we go 
through the argument backwards from Cx^/Cy-^ = F, to see that the triangle is also well defined. See 



1 38] for alternative discussion. 



Quadrilateral computations : Again we proceed in the spherical case and later adapt the hyper- 
bolic formulas. Consider a short quadrilateral with vertices and midpoints mi = mid(ai, Oj+i), 
s.t. distance{ai, Uj) < vr , Vi,j G {l,---4}. This means that not only the sides, but also the di- 
agonals are short. Now, denote by mo the midpoint of the diagonal (01,03) and by yo its half 
length. Similarly, denote by yi = ^distance{ai,ai+i). Then, each of the triangles {01,02,03) 
and (a3,a4,ai) are uniquely determined by their midpoints (mo, mi, 7712) and (mo,m3,m4), which 
shall be denoted by A12 and A34, respectively, with the same notation referring to their respective 
areas. Further, we denote Xij = distance{mi,mj). From the triangular analysis, we know that 
Cxi2 /Cyo = CA12/2, Cx3^/Cyo = Ca34/2- Similarly for the other partition, C^js/Cy^, = C'a23/2, 
Cx4i/Cyi^ = Ca4i/2i where y^ = ^ distance (02,04), m^ being its midpoint, and so on. There- 
fore, Cxj^2/Cx34 = C'ai2/2/Ca34/2 , Cx23/Cx4i = C ^23 /2 / C /2 ■ These equations generalize the 



CENTRAL VARIATIONAL PRINCIPLE 



27 



parallelogram relation on the plane, xi2/x^i = 1 = X23/a;4i, but contrary to the plane, they impose 
no constraint on the midpoints. To compensate for this fact, in the spherical case the 4 midpoints 
uniquely determine the area as well as the specific geometry of the short quadrilateral. In other words 
mo = mo({mi}) and similarly for mg. To see this, let's denote tq = ^Q^iai^a^). le. tq = (mo,uo), 
l^^ol = Vo- Also, we denote by TZm the involution through m and consider the element of 50(3) defined 
as cj^2 ■= ^mi^m2- Similarly, a^^ := TZm-^lZm^. Now, for any element a G 5*0(3), consider the vector 
field Xfj C TS"^ defined by r G X^j iff ^{t) = (m_,m+) s.t. m+ = cr^(m_). Then, the condition 
which guarantees that triangles A12 and A34 compose into a quadrilateral □1234 can be written as 
To G -X'o-12 n But, To G -^^0-12 onty if the pole pi2 of U12 lies in the polar line of ro, which is 

defined as the orthogonal geodesic to at "^o- Thus, the first condition obtained is that mo,pi2 and 
P34 be collinear. That is, Det[mo, Pi2: P34I = 0. In fact, we have more: pi2 = rhi x m2, P34 = 

X 7714. Therefore Det\rhn,mi x 1112, tfiR x = 0. Wc can thus write mn = z/\z\, where z = 
-212(^1 X "^2) + -234(^^3 X 7714). Further: denoting /ii2 =distance (^0,^12) and /X34 =distance (^0,^*34), 
we have {Sf^^^T^^^l = |Tj^J = IS'^g^T^gJ, but since {S^^^^^xiil = I'^o x (^1 x '^2)! and \S^^^Sj:^^\ = 
\rho X (7713 X 771,4)1, it follows that \rho x (mi x 7772) | • |m3 -7714! = |mo x (7773 x 7774)! " ""^2 1- Then, we 
get: 7770 = z/\z\, where z = aO^r^jl^s x 7774) + /3Ca;34 (7711 X 7772) , a,/3 G {—1, 1}. The sign choices must 
be made with care, but for small convex quadrilaterals we take the + choice twice. In this case, with 
R = C:ci2Cx3i + Cx23Cxii - Cxi3Cx24 5 ^6 have that \z\^ = = C^^^ + Cl^^ - 2C:E^^C:c3iR- Analogous 
equations hold for 777^. We have fixed the geometry of the small convex quadrilateral uniquely from 
the midpoints {mi, • • • , m4}, as mentioned earlier. The exception is when TZ\R.2'T^^T^A = 1 , but now, 
contrary to the euclidean plane, this is really the exception, not the rule. As for the area, from (6.8) we 
have S'ai2/2 = \ {Cxsa - RCx^} > ^'^34/2 = \ {Cxu - RC^^^}. Then, in this simpler case, denoting 
the quadrilateral area by P4, we get OP4/2 = R , P4 = 2o-i234 • C~^{Cxr2Cxa4 + CxisCx^ - Cx^Cx^a] ■ 
Here we have restricted to convex quadrilaterals decomposed in small triangles (area < vr). Greater 
care is needed for the sign choices in the expression of 777,0, otherwise. As for the restrictions, in this 
simpler case, by imposing ttiq • ttt., > 0, we get Z)i23 > 0, and cyclic . Again, in transposing to the 
hyperbolic plane, we don't have as many worries about sign choices, however in this case the analogous 
area function, P4 = 2a"i234 • C~^{R), where R = C12C34 + 023041 — O13O24, is well defined only when 
1^1 < 1 and, as in the triangular analysis, once this holds the convex quadrilateral exists. And so on 
for the general case. 
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